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INTRODUCTION 


This  project  was  initiated  by  Larry  Gedney  of  the  Geophysical 
Institute  of  the  University  of  Alaska.  About  halfway  through 
the  study  Larry  Gedney  retired,  and  David  Stone  and  John  Davies 
became  the  principal  investigators. 

The  original  goal  of  the  proposal  was  to  use  the  Dali  City 
earthquake  sequence  of  February  and  March,  1985  (fig.l)  to 
investigate  the  crustal  structure  of  north-central  Alaska.  The 
scope  of  the  proposal  was  later  modified  to  include  studies  of 
a  group  of  earthquakes  just  north  of  the  Alaska  Range,  the  Gold 
King  earthquakes  (fig.l)  and  later  still  to  include  data 
gathered  from  chemical  explosions  detonated  by  the 
U. S . Geological  Survey  as  part  of  their  Trans-Alaska  Crustal 
Transect  (TACT)  program  (fig. 3).  When  we  embarked  on  the  study 
of  the  Yukon  Flats  area  using  the  TACT  explosions,  there  were 
insufficient  funds  left  in  the  AFOSR-85-0266  contract,  so  we 
supplemented  these  funds  through  a  proposal  to  the  National 
Science  Foundation  (NSF)  and  through  the  use  of  student 
interns,  also  funded  by  the  NSF. 

As  with  any  project  of  this  size  and  scope,  it  overlapped  with 
other  ongoing  studies  of  the  crustal  structure  of  Alaska,  most 
notably  a  study  of  a  Fairbanks  -  Fort  Yukon  -  Barter  Island 
line,  and  a  study  of  Interior  Alaska  between  Fairbanks  and  the 
Alaska  Range.  The  basic  funding  for  these  latter  two  studies 
came  from  a  combination  of  State  of  Alaska  funds,  the 
U.S. Geological  Survey,  and  the  Rice-University  of  Alaska 
Industrial  Associates  Program,  with  considerable  indirect 
support  from  this  contract  ( AFOSR  85-0266)  in  terms  of  such 
things  as  the  analysis  of  common  data  sets. 


RESEARCH 

Although  all  the  work  completed  has  the  same  common  goal  of 
investigating  the  crustal  structure  of  Alaska,  it  divides  into 
roughly  five  parts. 

Dali  City  Events. 

The  occurrence  of  the  Dali  City  sequences  of  events,  which 
started  in  February  of  1985  (Estabrook  et  al,  1985;  Dixon, 

1989)  and  ended  (?)  with  a  magnitude  5.5  event  and  aftershocks 
in  July  of  1987,  allowed  the  crustal  models  being  used  for 
earthquake  locations  to  be  tested  and  refined.  A  description 
of  these  events  and  the  focal  mechanism  studies  can  be  found  in 
Appendix  I,  ("seismicity  and  focal  mechanism  studies  of  the 
1985  Dali  City  earthquakes,  north-central  Alaska"  by  Estabrook 
et  al.)  The  initial  crustal  models  deduced  from  these  events 
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are  described  by  Gedney  in  his  annual  technical  report  for  this 
contract  (Appendix  II,  Crustal  structure  studies  utilizing 
earthquake  sequences,  L. Gedney)  and  a  more  detailed  comparison 
made  with  other  crustal  models  in  Estabrook,  1985  (abstract 
given  as  Appendix  III),  and  Estabrook  et  al.  1988,  (Appendix 
IV)  . 

Gold  King  events. 

The  Gold  King  events,  so-called  because  of  their  proximity 
to  the  Gold  King  airstrip  and  seismometer  station  (fig.  1)  have 
considerable  importance  to  our  understanding  of  the  crustal 
structure  of  Alaska.  Their  importance  lies  in  the  fact  that 
they  apparently  mark  the  northeasternmost  seismic  activity 
directly  associated  with  the  downgoing  Pacific  plate  (Gedney 
and  Davies,  1986,  Appendix  V).  The  group  of  events  in  question 
lie  at  depths  between  100  and  140  km,  and  by  combining  them 
with  the  hypocenter  contours  for  events  to  the  west  and  south, 
the  location  of  that  portion  of  the  underthrusting  Pacific 
plate  can  be  determined.  To  the  east  of  the  Gold  King  events, 
it  appears  that  the  Pacific  plate  is  absent.  There  are  some 
hypocenters  about  300  km  south,  south  of  the  Wrangell 
Mountains,  that  are  consistent  with  the  Pacific  plate  reaching 
100km  depths  beneath  the  active  volcanic  center  of  Mt.  Wrangell 
itself,  thus  requiring  a  break  or  a  sharp  bend  in  the  Pacific 
plate  (Stephens  et  al.  1984?  Perez  and  Jacob,  1980).  If  it  is 
assumed  that  the  boundary  between  the  Pacific  plate  near  Gold 
King  and  where  it  is  observed  to  be  dipping  beneath  the 
Wrangell  Mtns  further  to  the  south  is  a  transform  boundary, 
then  this  boundary  can  be  depicted  by  a  small  circle  about  the 
Pacific  -  North  American  relative  motion  pole,  as  shown  in 
figure  2.  The  shaded  area  shown  in  figure  2  is  thus  the  area 
where  the  surface  geology  of  Alaska  is  underlain  by  oceanic 
crust.  The  interaction,  or  coupling  between  Alaska  and  the 
underthrust  Pacific  plate  may  be  generating  the  stresses 
causing  the  intraplate  seismicity  of  Interior  and  Arctic  Alaska 
as  well  as  producing  complex  crustal  velocities  in  southern 
Alaska . 

Yukon  Flats  Experiment 

The  Yukon  Flats  experiment  was  designed  to  xp'oit  the 
fact  that  the  U.S. Geological  Survey  TACT  program  .nated  a 
number  of  large  (up  to  6000  lb)  chemical  explosions.  The  basic 
aim  of  our  experiment  was  to  determine  the  crustal  structure 
beneath  the  Flats  (Technical  report  by  Stone  and  Davies,  1987, 
Appendix  VI) . 

The  geologic  and  tectonic  history  of  the  Yukon  Flats  is 
very  poorly  understood.  It  is  located  at  a  junction  of  several 
terranes  and  major  fault  systems  including  the  Tintina  and 
Kalt^g  faults  and  the  Porcupine  lineament  (figures  1,3,4). 
Because  there  are  few  constraints  on  its  history,  it  has  been 
postulated  to  be  underlain  by  ocean  crust  as  a  result  of  either 
entrapment  or  extension,  and  has  also  been  postulated  to  be 
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down-dropped  continental  crust.  The  basin  itself  is  known  to 
be  deep,  with  oil  industry  estimates  of  15  to  20  thousand  feet 
of  sediment  on  top  of  local  basement.  In  1987  the  TACT  program 
detonated  all  the  large  shots  shown  in  figure  5,  plus  a  number 
of  smaller  ones.  The  TACT  group  has  plans  for  the  early  1990s 
to  continue  the  shots  to  the  north  along  the  Trans-Alaska 
pipline  route  as  well  as  running  an  east-west  line  between 
Manley  and  Circle  hot  springs  (figure  1).  With  the  future 
shots,  we  hope  to  be  able  to  add  a  shot  off  the  east  end  of  the 
Yukon  river  stations  shown,  and  thus  more-or-less  reverse  the 
profile.  At  the  present  time  we  only  have  data  from  the  1987 
shots,  which  thus  gives  us  an  unreversed  profile  along  the 
Yukon  river,  but  also  gives  us  a  fan-like  distribution  of  paths 
across  the  Yukon  Flats  for  a  more  general  analysis  of  velocity 
versus  depth  (figure  5). 

For  the  1987  investigation  of  the  Yukon  Flats  area,  all  of 
the  stations  were  temporary  with  the  exception  of  Fort  Yukon, 
which  is  part  of  the  University  of  Alaska  permanent  network. 

Out  of  the  13  temporary  stations,  all  but  2  recorded  at  least 
some  of  the  shots.  The  best  records  were  those  that  were 
recorded  on  either  the  borrowed  USGS  5-day  recorders,  or  those 
that  transmitted  directly  to  Fairbanks  via  telephone. 

Figure  5  shows  the  travel-time  plots  for  the  Yukon  Flats 
fan  shooting  experiment.  The  interpretation  of  these  data  was 
greatly  aided  by  a  serendipidous  earthquake  at  Dali  City.  The 
magnitude  5.5  event  was  large  enough  for  us  to  identify  the  Pn 
arrivals  with  ease,  and  thus  to  "calibrate"  ourselves  to 
recognize  Pn  in  the  much  less  energetic  explosion  data.  From 
the  travel  times  we  deduce  that  the  crustal  structure  beneath 
the  Yukon  Flats  and  adjacent  Yukon-Tanana  uplands  consists  of 
an  8.7  km  thick  layer  with  a  velocity  of  5.4  km/s,  a  22.6  km 
layer  with  a  velocity  of  6.0  km/s,  and  a  Moho  velocity  of  8.45 
km/s.  This  gives  a  total  crustal  thickness  of  33.3  km,  which 
is  in  accord  with  the  models  derived  by  Estabrook  (1985)  and 
Estabrook  et  al.,  (1988). 

Permanent  network/Explosion  data. 

The  same  TACT  explosions  were  recorded  on  the  permanent 
network  operated  by  the  Geophysical  Institute  of  the  University 
of  Alaska  (figures  1  and  7).  To  investigate  the  crustal 
structure,  a  composite  travel-time  curve  has  been  assembled. 

The  curve  includes  data  along  a  line  from  near  to  the  Canadian 
border  on  the  Alcan  highway  to  the  Koyukuk  river  on  the  Trans- 
Alaska  pipline  route.  The  travel-time  curve  is  shown  in  figure 
8.  Although  all  of  the  composite  line  is  in  what  is  loosely 
termed  Interior  Alaska,  and  most  of  it  traverses  the  Yukon- 
Tanana  composite  terrane,  it  is  probable  that  it  samples 
several  different  crustal  configurations.  Taking  the  travel¬ 
time  curve  at  its  face  value,  it  gives  a  crustal  structure 
consisting  of  a  2.3  km  thick  layer  with  a  velocity  of  5.2  km/s, 
a  21.7  km  layer  with  a  velocity  of  5.9  km/s  and  a  Moho  velocity 
of  7.8  km/s.  This  gives  a  total  crustal  thickness  of  24  km, 


considerably  thinner  than  that  predicted  by  most  models 
(Estabrook  et  al.,  1988;  Stone  et  al . ,  1987). 

Another  approach  to  modeling  the  TACT  explosion  data  was 
taken  by  Johnson,  1989,  (abstract  given  as  Appendix  VII),  who 
treated  the  explosions  as  earthquakes,  and  located  them  using 
the  standard  HYPOELIPSE  reduction  program,  but  with  various 
crustal  models.  The  general  conclusions  from  her  study  are 
that  for  the  Fairbanks  area,  on  which  she  concentrated,  the 
location  parameters  were  not  consistently  better  for  any  one 
model  over  any  other. 

In  our  annual  report  dated  September  1987  (Appendix  VII) 
Stone  and  Davies  described  a  series  of  explosions  proposed  for 
the  Fort  Greely  area  that  we  hoped  to  use  for  further  crustal 
structure  experiments.  These  explosions  were  designed  to 
create  tank  traps,  and  the  combination  of  poor  coupling  (they 
were  designed  to  blow  dirt  out  rather  than  put  energy  into  the 
ground)  and  bad  weather  conditions  has  prohibited  us  from 
getting  useable  data  so  far. 

Barter  Island  -  Fairbanks  refraction 

Another  approach  to  large  scale  crustal  modeling  is  that 
of  Estabrook  (1985)  and  Estabrook  et  al.,  (1989).  In  this  case 
well  located  earthquakes  near  Barter  Island,  Fort  Yukon  and 
Fairbanks  were  used  as  energy  sources,  and  arrival  times  at  the 
seismometer  arrays  at  Barter  Island,  Fort  Yukon  and  Fairbanks 
used  to  construct  travel-time  curves.  As  can  be  seen  from 
figures  2  and  3  and  Table  1.  from  Appendix  IV,  the  depth  to 
basement  varies  from  about  33km  in  the  Fairbanks  area  through 
35  km  near  Dali  City  to  about  40  km  beneath  the  Brooks  Range, 
with  Moho  velocities  between  7.88  and  7.94  km/s. 


CONCLUSIONS 

Figure  9  shows  a  composite  of  crustal  models  for  Alaska.  It 
must  be  recognized  that  these  models  cover  most  of  Alaska,  and 
thus  large  differences  are  expected.  In  general  it  can  be  seen 
that  the  models  for  the  areas  north  of  Fairbanks  are  reasonably 
similar  which  is  consistent  with  the  fact  that  relocations  of 
events  using  these  different  models  do  not  produce  any  clear 
preference.  South  of  Fairbanks  the  situation  appears  to  be 
much  more  complex,  as  might  be  expected  with  the  greater 
complexity  of  young  accreted  terranes  combined  with  a  currently 
underthrusting  Pacific  plate. 
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FIGURES 


Figure  1. 

This  figure  shows  the  setting  of  the  Yukon  Flats  (northerly 
shaded  area)  and  the  rest  of  Interior  Alaska  (loosely  defined 
as  including  the  Tanana  valley  (shaded)  northward  to  the  Brooks 
Range  on  the  north  boundary  of  the  map) .  The  small  3-letter 
codes  represent  the  tectono-stratigraphic  terranes  of  Jones  et 
al.,  (1981);  Howell,  (1985);  and  Churkin  et  al.,  (1982).  Also 
shown  are  the  larger  of  the  TACT  shot  points  (*)  with  the  shot 
point  number  and  explosive  charge  in  1000s  of  pounds.  The 
seismic  stations  used  are  shown  as  triangles,  solid  for 
permanent  and  open  for  temporary.  The  road  system  is  shown 
with  a  dot  and  dash  symbol.  The  square  labelled  GKC  in  the 
Tanana  valley  is  the  Gold  King  site  used  to  locate  the  northern 
edge  of  the  Pacific  plate.  The  M  in  the  west  central  part  of 
the  map  represents  Manley  Hot  springs,  and  the  temporary 
station  CS  Circle.  These  two  locations  approximate  the  ends  of 
the  proposed  1990  or  1991  TACT  lines. 


Figure  2 . 

The  Pacific-Alaska  plate  interaction  in  the  southern  Alaska, 
redrawn  from  Stone,  (1983).  Hypocenter  depth  contours  (in 
kilometers)  and  the  shading  pattern  represent  the  extent  to 
which  the  Pacific  plate  has  been  underthust  beneath  Alaska. 
Stars  represent  recent  volcanoes,  and  the  thin  lines  the  major 
fault  systems.  The  RM1  lines  are  the  hypothetical  tracks  of 
the  Pacific  plate  about  the  Euler  poles  of  relative  motion 
determined  by  Minster  et  al.,  (1974). 


Figure  3 . 

The  route  of  the  TACT  transect  is  shown  by  the  heavy  dashed 
line.  The  major  faults  are  shown  as  solid  lines,  but  it  should 
also  be  noted  that  all  the  terranes  are  fault  bounded.  The 
Yukon  Flats  can  be  seen  to  lie  in  the  complex  area  between  the 
Tintina  and  Kaltag  faults,  the  Porcupine  lineament,  and  the 
composite  terranes  of  the  Brooks  Range. 

Figure  4 . 

This  figure  shows  the  relationship  between  the  Yukon  Flats  and 
the  various  geologic  provinces  of  Alaska.  The  part  of  Interior 
Alaska  investigated  by  the  studies  reported  here  include  parts 
of  the  Ruby  Schist  belt,  the  Yukon  Flats,  the  Yukon-Tanana 
crystaline  province  and  part  of  a  complex  of  Paleozoic  through 
Cenozoic  carbonates  and  clastic  rocks. 


Figure  5. 


An  outline  of  the  Yukon  Flats  area  (shaded)  with  representative 
paths  from  the  shot  points  to  seismometer  stations  (triangles) 
shown.  The  shot  points  r.re  labelled  with  their  number  and  the 
size  of  the  shot  in  1000s  of  lbs.  The  location  of  the 
fortuitous  1987  Dali  city  event  is  also  shown. 

Figure  6. 

Travel-time  plots  for  all  the  explosion  data  recorded  by  the 
stations  in  the  Yukon  Flats  area  together  with  the  arrivals 
from  the  Dali  City  event.  Because  of  the  shot  point-station 
distribution,  these  data  represent  a  fan-shooting  experiment. 

Figure  7 . 

The  distribution  of  all  the  TACT  shot  points  (triangles)  many 
of  which  were  used  for  several  seperate  explosions.  All  the 
stations  used  in  the  various  studies  are  represented  by  stars. 
The  solid  lines  represent  the  vertical  plane  onto  which  all  the 
adjacent  shot  points  and  receivers  were  projected  to  produce 
the  travel-time  curve  shown  in  figure  8. 

Figure  8. 

A  travel-time  curve  constructed  using  all  the  data  along  the 
solid  line  shown  in  figure  7.  The  best  estimates  of  the 
velocities  are  also  shown.  Since  this  is  a  composite  of  shot 
points  and  stations  distributed  along  the  whole  length  of  the 
line,  it  is  effectively  a  reversed  profile. 

Figure  9. 

A  composite  of  the  various  crustal  velocity  models  for  Alaska. 
The  two  new  models  associated  with  this  study  are  shown  by  the 
heavy  lines  labelled  "TACT  1987". 
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Figure  1 . 

This  figure  shows  the  setting  of  the  Yukon  Flats  (northerly 
shaded  area)  and  the  rest  of  Interior  Alaska  (loosely 
defined  as  including  the  Tanana  valley  (shaded)  northward  to 
the  Brooks  Range  on  the  north  boundary  of  the  map) .  The 
small  3 -letter  codes  represent  the  tectono-stratigraphic 
terranes  of  Jones  et  al.,  (1981);  Howell,  (1985);  and 
Churkin  et  al.,  (1982).  Also  shown  are  the  larger  of  the 
TACT  shot  points  (*)  with  the  shot  point  number  and 
explosive  charge  in  1000s  of  pounds.  The  seismic  stations 
used  are  shown  as  triangles,  solid  for  permanent  and  open 
for  temporary.  The  road  system  is  shown  with  a  dot  and  dash 
symbol.  The  square  labelled  GKC  in  the  Tanana  valley  is  the 
Gold  King  site  used  to  locate  the  northern  edge  of  the 
Pacific  plate.  The  M  in  the  west  central  part  of  the  map 
represents  Manley  Hot  springs,  and  the  temporary  station  CS 
Circle.  These  two  locations  approximate  the  ends  of  the 
proposed  1990  or  1991  TACT  lines. 


40* 


*50* 


Figure  2. 

The  Pacif ic-Alaska  plate  interaction  in  the  southern  Alaska, 
redrawn  from  Stone,  (1983).  Hypocenter  depth  contours  (in 
kilometers)  and  the  shading  pattern  represent  the  extent  to 
which  the  Pacific  plate  has  been  underthust  beneath  Alaska. 
Stars  represent  recent  volcanoes  and  the  thin  lines  the 
major  fault  systems.  The  RM1  lines  are  the  hypothetical 
tracks  of  the  Pacific  plate  about  the  Euler  poles  of 
relative  motion  determined  by  Minster  et  al.,  (1974). 
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bounded.  The  Yu  con  1 xats  can  be  seen  to  lie  in  the  complex 
area  between  the  Xintina  and  Kaltag  faults,  the  Porcupine 
lineament,  and  the  composite  terranes  of  the  Brooks  Range. 
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Figure  4. 

This  figure  shows  the  relationship  between  the  Yukon  Flats 
and  the  various  geologic  provinces  of  Alaska.  The  part  of 
Interior  Alaska  investigated  by  the  studies  reported  here 
include  parts  of  the  Ruby  Schist  belt,  the  Yukon  Flats,  the 
Yukon-Tanana  crystaline  province  and  part  of  a  complex  of 
Paleozoic  through  Cenozoic  carbonates  and  clastic  rocks. 


Figure  5. 

An  outline  of  the  Yukon  Flats  area  (shaded)  with 
representative  paths  from  the  shot  points  to  seismometer 
stations  (triangles)  shown.  The  shot  points  are  labelled 
with  their  number  and  the  size  of  the  shot  in  1000s  of  lbs. 
The  location  of  the  fortuitous  1987  Dali  city  event  is  also 
shown. 
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APPENDIX  I 


Seismicity  and  Focal  Mechanism  Studies 
of  the  1985  Dali  City  Earthquakes, 

No r th - cent r al  Alaska 


By  Charles  H.  Estabrook,  John  N.  Davies, 
Larry  Gedney,  Steven  A.  Estes  and  James  P.  Dixon 
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Abs tract 

In  February  and  March  of  1985,  several  moderate 
earthquakes  occurred  in  central  Alaska,  42  km  north  of  the 
Yukon  River  crossing  of  the  Trans-Alaska  Pipeline.  The  largest 
of  these,  on  9  March  (M  6.0)  was  felt  widely  in  Interior 
Alaska;  38  events  had  magnitudes  greater  than  4.0.  Several 
short  period  vertical  seismometers  were  installed  in  the 
epicentral  region  by  personnel  from  the  Geophysical  Institute 
of  the  University  of  Alaska.  A  P-  and  S-wave  seismic  velocity 
model  determined  using  several  wel 1  -  recorded  aftershocks 
indicates  that  granite  is  the  major  component  of  the  upper 
crust  in  the  epicentral  area.  This  velocity  model  was  used  to 
locate  the  aftershocks.  The  aftershocks  fall  into  two 
intersecting  trends.  The  main  trend  (strike  28°)  is  25  km  long 
and  dips  75-80°  E,  with  depths  ranging  from  the  surface  to  20 
km.  The  secondary  trend  strikes  43°.  These  aftershocks  are 
north  of,  but  spatially  separated  from,  the  aftershock  zone  of 
the  1968,  Ms  6.5  Rampart  earthquake.  They  spread  to  the  north 
and  south  along  assumed  fault  planes  related  to  the  largest 
events  and  along  a  possible  intersecting  fault,  with  the 
largest  earthquakes  in  the  sequence  occurring  near  the 
intersection  of  the  two  trends.  Focal  mechanism  solutions  of 
two  of  the  largest  events  show  left-lateral  strike-slip 
faulting  on  NNE  trending  planes. 
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Introduction 

In  February  and  March  of  1985,  a  sequence  of  five  moderate 
earthquakes  (magnitude  5  to  6)  occurred  in  central  Interior 
Alaska,  approximately  40  km  north  of  the  Yukon  River  crossing 
of  the  Trans-Alaska  Pipeline  (Fig.  1).  On  14  February  (0504 
UT),  an  event  of  magnitude  5.4  ( M^ )  occurred  which  was  followed 
by  a  normal  aftershock  sequence.  Figure  2  shows  that  activity 
had  nearly  ceased  by  the  beginning  of  March.  A  remote  seismic 
station  (BEH,  Fig.  3),  installed  on  17  February,  was  to  be 
removed  at  about  the  time  a  second  ma inshock - afte r shock 
sequence  began.  Early  in  the  morning  of  9  March  (4:58  am  AST; 
1358  UT)  an  event  of  magnitude  5.0  ( ML )  awoke  a  few  people  in 
the  Fairbanks  area.  This  was  followed  at  1408  UT  by  a 
magnitude  6.0  (ML)  shock,  the  largest  earthquake  to  occur  in 
Interior  Alaska  since  the  1968  Rampart  earthquake  (Gedney  et 
al . ,  1969b;  Huang  and  Biswas,  1983). 

The  1408  event  was  felt  widely  in  Interior  Alaska  --  as 
far  away  as  Delta  Junction,  300  km  southeast,  and  Coldfoot,  200 
km  north  in  the  Brooks  Range.  In  Fairbanks,  180  km  southeast 
of  the  epicenter,  the  Intensity  of  the  earthquake  ( Modified 
Mercalli  Scale)  was  IV;  in  Manley  Hot  Springs,  90  km  from  the 
epicenter,  where  grandfather  clocks  stopped,  the  Intensity  was 
V;  and  at  the  Yukon  River  Crossing,  where  an  Alaska  Department 
of  Transportation  Trailer  was  shaken  off  its  supports,  the 
Intensity  was  probably  VI. 


-4- 


In  the  epicentral  area,  where  there  are  no  manmade 
structures,  the  effects  of  the  quakes  included  numerous  stream 
overflows  and  spring  rejuvenation.  However,  there  was  no 
observed  ground  breakage.  Muddy  water  appeared  to  have 
overflowed  on  the  West  Fork  of  the  Dali  River  approximately  10 
km  due  north  of  the  epicentral  area  (Fig.  3).  Recent  stream 
icings  (spring  rejuvenation  following  the  largest  events)  were 
seen  along  the  north-south  portion  of  the  Dali  River  (F.  Weber, 
pers.  comm.)  and  along  the  east-west  portion,  west  of  Dali  City 
(the  site  of  a  prospectors'  camp  abandoned  in  1901).  These 
effects  were  observed  during  several  flights  over  the 
epicentral  area  by  helicopter  and  fixed-wing  aircraft  from 
February  through  April. 


Geologic  Setting 

The  Dali  City  earthquakes  occurred  in  a  region  of  low 
hills  between  the  Ruby  Mountains  and  the  Yukon  Flats  in  north- 
central  Alaska  (Fig.  1).  Pelitic  schists  and  granites  of 
Paleozoic  and  Mesozoic  age  characterize  the  geology  of  the 
epicentral  area  (Brosge  et  al . ,  1973).  Though  structural 
features  are  often  obscured  by  the  extreme  weathering  of  the 
schists  and  granites,  several  faults  have  been  mapped  in  the 
epicentral  area.  The  Dali  Mountain  fault  (Brosge  et  al . ,  1973; 
Brogan  et  al . ,  1975;  and  Plafker,  unpubl .  data)  trends  north- 
south  approximately  10  km  east  of  the  epicenters.  Several 
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northeast  trending  faults  have  been  mapped  in  the  area  to  the 
west  of  the  epicenters  (Dover  and  Miyaoka,  1985).  The  Tintina 
and  Kaltag  fault  systems  are  found  to  the  southeast  and 
southwest  of  the  epicenters  respectively.  These  faults  exhibit 
recent  movement  (Brogan  et  al .  ,  1975  ;  Estabrook  et  al .  ,  1985). 
An  aeromagnetic  lineation,  the  Porcupine  lineament  (Brosge  et 
al  .  ,  1970),  crosses  the  area  from  the  northeast.  To  the  east 
of  the  epicenters  is  the  Yukon  Flats,  a  Quaternary  basin  which 
may  have  developed  as  an  extentional  graben  from  strike-slip 
movement  along  the  Tintina  fault  (Kirschner  et  al . ,  1985). 

Data  Sources 

The  seismic  data  for  the  Dali  City  events  recorded  by  the 
Geophysical  Institute  can  be  divided  into  five  sets  (Table  1) 
based  on  differing  combinations  of  local  recording  stations 
used  to  supplement  the  existing  regional  network  (Fig.  1). 

Prior  to  the  installation  of  a  short-period  vertical 
seismometer  at  the  ALASCOM  Microwave  Repeater  Site,  Bench 
(BEH),  on  February  17,  the  closest  station  was  at  Indian 
Mountain  (IMA),  about  170  km  west  in  the  Koyukuk  Basin.  Hence 
the  first  data  set  was  recorded  with  no  local  stations.  BEH 
was  the  only  local  station  operating  until  two  additional 
stations  began  operating  on  March  15.  Data  recorded  with  the 
single  local  station,  BEH  and  the  regional  network  constitutes 
the  second  data  set.  Attempts  were  made  to  install  a  station 
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at  Fish  Microwave  Repeater  Site  (FISH,  Fig.  3)  in  February  but 
the  effort  was  abandoned  due  to  extreme  winter  weather 
conditions  (extremes  greater  than  -41°C  and  75  km/hr  winds  were 
recorded  at  BEH  the  same  day)  and  a  decrease  in  aftershock 
activity  (Fig.  2).  By  15  March  attempts  to  install  this 
station  were  successful  and  data  were  being  collected  from  BEH, 
FISH  and  a  third  local  station  emplaced  at  a  remote  gate  valve 
on  the  Trans-Alaska  Pipeline  (RGV,  Fig.  3).  Data  from  those 
stations  in  combination  with  that  from  the  regional  network 
comprise  the  third  data  set.  The  local  network  was  modified 
slightly  in  April  (fourth  set)  and  entirely  removed  in  June 
1985  (fifth  set ) . 

Data  from  all  of  the  Interior  Alaska  stations  used  by  the 
Geophysical  Institute  are  recorded  in  analog  form  on  a  Teledyne 
Geotech  Deve loco rde r .  In  addition,  an  event  detection  system, 
developed  by  Steven  Estes,  records  signals  from  the  same 
stations  in  digital  form  (8  bits/sample  at  50  samples/s).  An 
earthquake  is  recorded  when  the  ratio  of  a  short  term  average 
of  the  signal  to  the  long  term  average  exceeds  3.5  at  three 
stations  during  a  thirty  second  period.  Both  analog  and 
digital  data  were  used  in  this  study. 

Focal  Mechanisms 

Fault  plane  solutions  were  determined  for  the  two  largest 
earthquakes  (14  February  0504  and  9  March  1408)  in  the 
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sequence.  In  addition,  a  solution  was  obtained  for  an  earlier 
earthquake  in  the  same  area  which  occurred  nearly  ten  years  to 
the  minute  before  the  largest  one  in  1985  (1419,  9  March, 
1975).  The  mechanisms  were  determined  using  a  linearly 
increasing  velocity  model  with  the  hypocentral  depths  fixed  at 
10  km.  The  parameters  of  these  focal  mechanisms  are  given  in 
Table  2.  All  three  events  are  strike-slip  solutions  with  very 
similar  nodal  plane  orientations  (Fig.  4).  The  14  February 
mechanism,  which  utilized  only  regional  Alaska  stations,  is 
very  dependent  upon  an  emergent  P-wave  arrival  at  Barrow  (BRW, 
the  dilatation  in  the  northwest  quadrant).  The  9  March 
solution  utilizes  first  motion  readings  from  the  Canadian 
Seismograph  Network,  the  U.  S.  Geological  Survey's  Central 
California  Network  and  the  Lamont -Doherty  Geological 
Observatory  station  at  Palisades,  NY.  This  mechanism  is  much 
more  tightly  constrained.  Note  that  a  high-angle  NNW  trending 
nodal  plane  is  common  to  both  solutions.  The  9  March,  1975 
event  occurred  25  km  to  the  south  of  the  other  two  and  was 
smaller  in  magnitude  (M  4.6)  (Estabrook,  1985). 

L 


Crustal  Velocity  Determination 

The  local  P-  and  S-wave  velocity  structure  in  the  Dali 

City  area  was  studied  utilizing  the  availability  of  digital 

data  and  the  close  station  spacing.  The  approach  taken  was  to 

determine  a  V  /V  ratio  for  the  Dali  City  area  and  then  use 
p  s  J 
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this  ratio  in  the  HYPOELLIPSE  earthquake  location  program 
(Lahr,  1982)  to  find  a  minimum  in  the  location  and  travel-time 
error  estimates  as  a  function  of  the  model  velocity. 

Figure  5  shows  an  example  of  a  Dali  City  event,  with  well 
recorded  P-  and  S-waves,  which  was  used  in  the  velocity 
studies.  Wadati  diagrams  were  used  to  determine  the 

Vp/Vs 

ratio  (Fig.  6).  The  resultant  V  /V  ratio  of  1.68  +  0.03  (for 

P  s 

nine  earthquakes)  is  very  close  to  1.69  that  was  determined  by 
Gedney  et  al  .  (1969a)  from  aftershocks  of  the  1967  Fairbanks 

earthquakes . 

Studies  of  seismic  velocities  in  granitic  and  metamorphic 

rocks  in  the  Ukrainian  Shield  (Lebedev  and  Korchin,  1979)  show 

that  Vp/Vs  is  higher  for  mafic  rocks  (1.77  +  0.06)  than  for 

silicic  rocks  (1.68  +  0.09).  Pelitic  schists  have  values 

comparable  to  silica-rich  igneous  rocks  (Press,  1966).  The 

presence  of  silica-rich  metamorphics  and  granites  at  the 

surface  in  the  Dali  City  area  (Brosge  et  al .  ,  1973),  combined 

with  a  V  /V  ratio  of  1.68,  indicates  that  quartz-rich  rocks 
P  ® 

(granites  or  schists)  occur  at  depth  in  the  region. 

In  an  attempt  to  determine  the  P-wave  velocity  for  the 
Dali  City  region,  the  crustal  P-wave  velocity  in  our  model 
(Table  3)  was  increased  at  0.1  km  increments  from  5.5  to  7.0 
km/s.  These  velocities  were  used  in  the  HYPOELLIPSE  program  to 
locate  10  earthquakes  for  which  very  good  P-  and  S-wave 
arrivals  were  recorded  at  the  local  stations  BEH,  RGV  and  FISH. 
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The  results  of  this  study  (Fig.  7)  show  that  as  we  increase  the 

2 

crustal  velocity  in  our  model,  the  location  errors  ( ( ERH  + 

2  1/2 

ERZ  )  )  and  root  -  mean - squa re  (RMS)  of  the  travel-time 

residuals  decrease.  It  was  hoped  that  a  minimum  in  location 
error  and  RMS  could  be  found.  However,  because  the  local 
station  distribution  is  nearly  linear,  there  is  a  trade-off 
between  epicentral  location  and  depth  as  the  velocity 
increases.  A  crustal  velocity  of  6.0  km/s  corresponds  to  the 
point  at  which  RMS  and  location  errors  asymptotically  approach 
a  lower  limit.  Though  velocities  higher  than  this  also  have 
very  small  location  errors  and  RMS,  velocities  greater  than  6.3 
km/s  yield  depths  that  are  unreasonable  for  Interior  Alaska 
(Fairbanks  area:  Gedney  et  al . ,  1980)  and  south  -  central  Alaska 
(Castle  Mountain  fault:  J.  Lahr,  pers.  comm.,  1984). 

Taking  into  consideration  the  factors  of  surface  geology, 
a  minimum  value  for  the  upper  crustal  velocity  and  the 

Vp/Vs 

ratio,  we  can  constrain  the  geology  and  velocity  at  depth  by 
making  the  assumption  that  if  earthquakes  in  the  Dali  City  area 
have  depths  similar  to  those  of  other  earthquake  sequences  in 
continental  Alaska,  then  the  velocity  of  the  upper  crust  must 
lie  between  6.0  and  6.3  km/s.  Since  the  stations  used  in  the 
local  velocity  study  vary  from  15  to  40  km  from  the  epicenters, 
the  rays  leaving  an  earthquake  at  10-15  km  depth  will  generally 
travel  upward  to  the  stations.  Thus  the  local  velocity 
corresponds  to  that  of  an  upper  crust  having  a  thickness  of  at 
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least  the  depths  of  the  earthquakes.  This  value  of  the  upper 
crustal  velocity  agrees  with  the  6.04  km/s  velocity  determined 
from  a  regional  earthquake  in  the  eastern  Brooks  Range 
(Estabrook  et  al . ,  1985). 

Press  (1966)  reports  that  for  depths  to  15  km  (5  kbars), 
granites  have  P-wave  velocities  of  about  6.3  km/s  while  schists 
have  velocities  that  vary  from  6.5  to  6.9  km/s.  Since  the 
upper  crustal  velocity  is  6.0-6. 3  km/s,  the  rock  type  of  the 
material  at  depth  (about  10  km)  must  then  be  mainly  granitic 
with  possibly  some  schist,  although  some  sedimentary  rock  types 
(graywacke  and  some  limestones)  also  fit  the  P-wave  velocity 
criteria  of  6. 0-6. 3  km/s.  More  mafic  rocks  (gabbro,  dunite  and 
ophiolitic  rocks)  are  ruled  out  as  major  components  of  the 
upper  crust  because  their  P-wave  velocities  are  too  high. 

Seismicity 

Because  four  distinct  station  distributions  were  used 
during  the  two  month  period  to  record  the  Dali  City 
earthquakes,  a  uniform  interpretation  of  the  sequence  is 
difficult.  The  first  data  set  (Fig.  8a)  contains  earthquakes 
located  only  with  stations  at  Pn  distances.  During  this  time, 
the  depths  of  the  events  are  not  well  constrained.  Since  BRW 
was  the  only  station  located  to  the  north  of  the  epicenters 
(further  than  600  km),  the  epicenters  will  more  likely  be 
mislocated  in  the  north-south  direction.  Lee  and  Stewart 
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(1981)  showed  that  when  there  is  a  large  gap  in  station 
distribution,  the  epicenter  will  be  mislocated  away  from  its 
real  location  toward  that  gap.  The  second  data  set  (Fig.  8b) 
contains  earthquakes  located  with  the  regional  stations  (as  in 
set  1)  and  a  single  local  station  to  the  south  of  the  Dali  City 
area  (BEH).  In  this  set,  errors  in  depth  are  reduced  because 
the  ratio  of  depth  of  the  event  to  the  distance  to  the  local 
station  is  about  1:2.  Epicentral  errors  are  still  likely 
because  of  the  gap  problem,  though  the  errors  should  be  smaller 
than  without  the  local  station.  Earthquakes  in  the  third  set 
(Fig.  8c)  are  located  with  the  regional  network  (same  as  the 
first  two)  and  three  local  stations.  With  the  linear  local 
station  distribution,  computed  hypocenters  will  move  toward  the 
largest  gap,  which  is  away  from  the  axis  of  the  local  stations. 
This  effect  is  minimized  with  accurate  S-wave  readings  for  the 
smaller  events  and  the  inclusion  of  regional  stations  FYU  and 
BM3  (Fig.  3)  to  the  east  of  the  epicenters  in  the  hypocentral 
location  process. 

Figure  8  contains  maps  of  all  the  solutions  obtained  with 
data  to  the  end  of  April  1985.  Notice  that  scatter  is  reduced 
when  local  stations  are  used  in  the  solutions.  The  trend  of 
epicenters  from  the  data  sets  delineates  a  N28°E  trend.  Figure 
8d,  in  addition,  shows  a  N43°E  trend  that  intersects  the  main 
trend  in  the  vicinity  of  the  largest  events. 
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To  delineate  the  fine  structure  of  the  aftershock  trend, 
the  earthquakes  located  with  three  local  stations  plus  the 
regional  network  were  graded  on  the  basis  of  quality.  It  is 
assumed  that  mislocated  earthquakes  will  have  large  location 
and  travel-time  errors.  From  Figure  9a,  which  shows  the  best 
solutions  of  this  period  (  ERH  and  ERZ  _<  10  km,  RMS  _<  0.5  s),  a 
very  clear  NNE-SSW  lineation  along  with  a  possible  ENE-WSW 
lineation  emerges.  Cross  -  sections  of  this  data  set  are  shown 
in  Figures  9b  and  9c.  Figure  9b  shows  that  the  trend  dips 
steeply  to  the  east.  Drawn  on  this  cross-section  is  the  north 
trending  nodal  plane  of  the  9  March  event.  Notice  that  the 
dips  are  nearly  the  same.  The  seismic  width  of  the  aftershock 
zone,  about  10  km,  is  likely  to  be  thinner  but  beyond  the 
resolution  possible  with  the  network  configuration.  in  the 
longitudinal  section  (Fig.  9c),  the  aftershock  zone  is  roughly 
20  km  long.  It  appears  in  plots  of  aftershocks  of  the  14 
February  event  prior  to  the  installation  of  BEH  (Fig.  8a). 

With  the  exception  of  the  Dali  Mountain  fault  about  10  km 
to  the  east,  there  is  no  structural  feature  in  the  epicentral 
area  with  topographic  expression  subparallel  to  the  N28°E  trend 
of  the  aftershocks.  The  N43°E  trend  is  parallel  to  an 
escarpment  extending  from  the  Ray  Mountains  in  the  west,  to  the 
east  along  the  northern  boundary  of  the  Yukon  Flats.  Also  on 
this  trend  are  several  hot  springs  in  the  Ray  Mountains,  which 
may  indicate  that  the  lineament  extends  through  the  granite  of 
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the  mountains.  On  a  geologic  map  of  the  Beaver  quadrangle 
(1  :  250000  :  Brosge  et  al .  ,  1973),  there  is  no  mapped  fault 
parallel  to,  or  geologic  di scont inui ty  across,  the  escarpment. 
However  recent  mapping  just  to  the  west  of  the  epicentral 
region  indicates  that  a  northeast  trending  fault  parallels  the 
escarpment  (Dover  and  Miyaoka,  1985).  This  fault  (with  an 
associated  1-3  m  wide  zone  of  fault  gouge)  was  observed  in  a 
quarry  off  the  Dalton  Highway  by  one  of  the  authors 
( Estabrook ) . 


Discussion 

Focal  mechanisms  of  the  two  largest  events  show  either 
left-lateral  faulting  on  NNE  oriented  planes  or  right  -  lateral 
faulting  on  ESE  oriented  planes.  The  major  trend  of  the 
aftershocks  is  parallel  to  the  NNE  trending  nodal  plane.  On 
this  basis,  and  in  consideration  of  the  parallelism  of  the  NNE 
nodal  planes-with  the  mapped  trace  of  the  Dali  Mountain  fault, 
the  NNE  plane  is  preferred. 

Figure  10  is  a  plot  of  18  of  the  largest  events  of  the 
sequence  (see  Table  4  for  location  parameters).  These  events 
were  selected  because  they  were  all  large  enough  to  be  recorded 
at  BRW  (Barrow),  a  low-gain  station  that  is  only  able  to  record 
events  with  magnitudes  greater  than  4.0  at  this  distance. 

Notice  that  the  events  cluster  near  the  intersection  of  the  two 
aftershock  trends  (compare  with  Figure  9a).  The  large  events 
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also  show  temporal  grouping  and  possible  westward  migration. 

The  events  of  14  February  and  9  March  cluster  in  the  east, 
those  of  10  March  are  located  3-5  km  to  the  southwest.  Either 
this  ic  a  migrat 1  on  on  a  single  fault  plane  oriented  northeast 
or  faulting  on  a  series  of  parallel  north-south  trending  faults 
which  intersect  the  secondary  aftershock  trend.  If  the  NNE 
nodal  plane  is  the  fault  plane  and  the  larger  shocks  actually 
do  migrate  from  east  to  west,  then  we  interpret  that  there  is  a 
parallel  set  of  north-south  trending  faults  in  the  vicinity  of 
the  Dali  Mountain  fault  that  intersect  the  northeast  trending 
faults . 

The  location  of  the  larger  earthquakes  in  a  relatively 
small  area  of  the  aftershock  zone  implies  that  stress  must  be 
concentrated  in  this  region.  The  probable  asperity  (stress 
concentration)  is  at  the  intersection  of  the  two  aftershock 
trends.  The  majority  of  the  aftershocks  occurred  on  a  trend 
parallel  to  the  fault  plane  of  the  9  March  mainshock.  The 
N43°E  aftershock  trend  would  then  represent  sympathetic  motion 
on  an  intersecting  fault. 

The  seismicity  data  compiled  by  Estabrook  et  al.  (1985) 
indicates  that  the  Dali  City  aftershocks  plot  just  north  of  the 
aftershock  zone  of  the  1968,  M^  6.5  Rampart  earthquake 
northward  to  66.5°N  with  slightly  different  trends  (Fig.  11). 
This  map  clearly  shows  that  since  the  two  earthquake  sequences 
do  not  touch  or  follow  the  same  trend,  they  occurred  on 
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different  lineaments.  The  Dali  City  aftershocks  are  parallel 
to  faults  to  the  west  (Dover  and  Miyaoka,  1985),  while  the 
Rampart  events  follow  a  trend  that  connects  two  segments  of  the 
Kaltag  fault  system  (Huang  and  Biswas,  1963),  stopping  just 
south  of  the  Yukon  River.  The  two  aftershock  zones  are 
separated  by  faults  of  the  Kaltag  fault  system. 

Portions  of  the  Kaltag,  Porcupine  and  Tintina  fault 
systems  are  all  found  within  the  region.  Both  the  Kaltag  fault 
and  the  Porcupine  lineament,  inferred  from  aeromagnetic  data, 
parallel  the  northeast  trending  escarpment.  To  add  to  the 
confusion,  a  line  can  be  drawn  which  follows  the  trace  of  the 
Tintina  fault  through  the  Dali  City  epicentral  region;  however, 
this  trend  is  not  consistent  with  either  of  the  two  aftershock 
trends.  Only  the  faults  mapped  by  Dover  and  Miyaoka  (1985)  to 
the  west  and  the  northeast  trending  escarpment  have  similar 
trends  to  the  aftershocks.  In  fact,  several  faults  that  they 
mapped  intersect  at  approximately  the  same  angle  as  that  of  the 
two  aftershock  trends  (15°).  On  this  basis  we  surmise  that  the 
Dali  city  earthquakes  occurred  at  the  intersection  of  two 
faults  of  the  same  system  as  those  mapped  by  Dover  and  Miyaoka, 
with  the  largest  earthquakes  occurring  near  the  intersection. 
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Figure  Captions 

Fig.  1.  (a)  Location  map  for  Dali  City  earthquake 

studies.  (b)  Interior  Alaska  seismograph  stations  (three 
letter  codes)  used  to  locate  the  Dali  City  earthquakes  of 
February  to  April  1985.  The  shaded  area  delineates  the 
location  of  all  of  the  seismicity  plots.  P5  and  P6  are  pump 
stations  along  the  Trans-Alaska  Pipeline.  The  Kaltag  and 
Tintina  faults  (heavy  lines  in  the  west  and  east,  respectively) 
and  the  Porcupine  lineament  (dotted  line)  are  from  King  (1969). 

Fig.  2.  Histogram  of  Dali  City  earthquakes  as  recorded  at 
Gilmore  Dome  (GLM). 

Fig.  3.  Map  of  surface  effects  of  the  March  1985 
earthquakes.  Hatchered  area  is  zone  of  spring  rejuvenation. 
Tick  marks  on  streams  indicate  zone  of  stream  overflow  and 
broken  ice.  The  asterisk  shows  the  location  of  the  M  6.0 

L 

event  of  1408,  9  March.  Solid  dots  are  earthquakes  with  M  > 
5.0.  The  Dali  Mountain  fault  is  from  Brosge,  et  al.  (1973). 

Fig.  4.  Fault  plane  solutions  of  the  14  February,  9  March 
1985  and  9  March  1975  earthquakes.  They  are  equal  area  lower 
hemisphere  projections  showing  quadrants  of  compressional  (dark 
circles)  and  tensional  (open  circles)  first  motion.  The 
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letters  'P'  and  'T'  are  the  pressure  and  tensional  axes 
respectively.  Tick  marks  on  the  nodal  planes  represent  the 
poles  to  the  other  planes. 

Fig.  5.  Seismogram  of  an  event  used  in  the  Vp/Vs  ratio  and 
crustal  structure  study. 

Fig.  6.  Wadati  diagram  of  the  earthquake  in  figure  5.  The 
line  connecting  the  points  is  a  least  -  squares  fit  giving  a 
Vp/Vs  ratio  of  1.68. 

2  2  1/2 

Fig.  7.  Location  error  ( ERH  +  ERZ  )  ,  travel-time  residual 

(RMS)  and  depth  as  a  function  of  P-wave  velocity  for  Vp/Vs  = 

1.68. 

Fig.  8.  Earthquakes  which  occurred  from  (a)  1  to  17  February 
1985,  prior  to  the  installation  of  BEH;  (b)  17  February  to  9 
March,  prior  to  the  1408  event;  (c)  9  March  to  the  installation 
of  FISH  and  RGV  on  15  March;  and  (d)  after  the  installation  of 
FISH  and  RGV  to  30  April.  Depths  are  indicated  by  letter:  A=0 
to  5  km,  B=5  to  10  km.... 

Fig.  9.  (a)  Earthquakes  recorded  during  the  period  from  15 

March  to  30  April  (as  in  Figure  8d)  with  ERH  and  ERZ  £  10  km, 
RMS  <  0.5  sec,  GAP  <  135°  and  no.  stations  >  5.  (b)  is  the  NW- 
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SE  cross  -  section  (A-A'J  of  the  Dali  City  earthquakes.  (c)  is 
the  NE-SW  cross-section  ( B  -  B '  )  . 

Fig.  10.  The  largest  events  (M  .  <  4.0)  in  the  Dali  City 
earthquake  sequence.  Their  location  indicates  the  location  of 
the  probable  asperity. 

Fig.  11.  Seismicity  of  nor th - cent ral  Alaska.  The  earthquakes 
plotted  here  are  from  the  catalog  compiled  by  Estabrook  et  al . 
(1985)  for  the  years  1967  to  1985.  They  have  ERH  and  ERZ  <  10 
km,  RMS  <  0.75  s  and  no.  stations  5.  The  earthquakes  before 
February  1985  are  plotted  as  pluses,  the  events  of  February- 
April  are  plotted  as  stop  signs.  Faults  are  from  Beikman 
(1980),  Dover  and  Miyaoka  (1985)  and  King  (1969).  The  number 
of  events  in  the  Rampart  aftershock  zone  is  not  representative 
of  the  actual  number  recorded  by  Alaskan  seismic  stations;  its 
extent,  however,  is  consistent  with  that  determined  by  Huang 
and  Bi swas  (  198  3)  . 
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Summary 


During  the  course  of  this  investigation,  two  earthquake 
sequences,  one  100  km  southwest  of  Fairbanks  and  one  180  km  to  the 
northwest,  were  studied  in  depth  to  enhance  the  present  knowledge  of 
crustal  structure  and  tectonic  plate  interactions  in  central  Alaska. 
The  former  sequence  consisted  of  intermediate-depth  earthquakes  at  the 
extreme  northern  tip  of  the  encroaching  Pacific  plate,  and  the  latter 
consisted  of  shallow  crustal  earthquakes  in  an  intra-plate  regime. 

The  deeper  cluster  to  the  south  was  found  to  be  a  product  of  down-dip 
tension  in  the  subducting  plate.  The  crustal  earthquakes  to  the  north 
result  from  left-lateral  faulting  which  arises  from  nearly  horizontal 
compressive  stress  from  the  south-southeast.  An  improved  crustal 
model  was  obtained  for  the  area.  The  best-fit  parameters  for  the 
central  part  of  the  zone  include  an  upper  mantle  P-wave  velocity  of 
7.91  km/sec  at  a  depth  of  34  km.  Overlying  this  there  are  two  crustal 
layers,  the  uppermost  of  which  is  24  km  thick  and  has  a  P-wave 
velocity  of  6.0  km/sec.  The  lower  crustal  layer  is  10  km  thick,  with 
a  velocity  of  7.2  km/sec.  Upper  mantle  (Pn)  velocities  were  observed 
to  vary  with  azimuth  around  the  margins  of  the  central  area  from  as 
low  as  7-63  km/sec  to  as  high  as  8.12  km/sec.  This  no  doubt  reflects 
the  flexing  and  deepening  of  the  Moho  beneath  the  root  of  the  Alaska 
Range . 


STATEMENT  OF  WORK 

One  objective  of  this  project  was  to  provide  additional 
information  on  the  stress  regime  in  a  very  complicated  area  of  Central 
Alaska  where  the  North  Pacific  plate  is  indenting  the  continent  in  the 

area  of  the  great  bend  of  the  Alaska  Range  around  Mt .  Mckinley.  This 

was  enhanced  by  the  discovery  of  a  tightly-compacted  cluster  of 
intermediate-depth  earthquakes  that  have  apparently  been  going  on 
unnoticed  for  an  indeterminate  period  of  time.  As  related  in  the 
following  section,  this  cluster  was  found  to  define  the  northernmost 
tip  of  the  encroaching  Pacific  plate. 

A  second  objective  was  to  obtain  an  improved  crust  and  upper 
mantle  model  for  Central  Alaska  in  a  region  about  which  little  was 

known.  This  was  made  possible  by  the  occurrence,  in  February,  March 

and  April  of  1985.  of  several  moderate  earthquakes  and  over  1,400 
aftershocks  in  an  area  previously  considered  to  be  relatively 
aseismic.  A  magnitude  6.0  mainshock  on  March  9  is  the  largest 
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earthquake  ever  to  have  been  documented  this  far  north  in  Alaska. 

Because  the  epicentral  region  was  in  close  proximity  to  the 
transportation  corridor  paralleling  the  Alaska  pipeline,  it  was 
possible  to  install  three  seismographic  stations  in  the  immediate 
vicinity,  and  to  telemeter  the  signals  directly  back  to  the 
Geophysical  Institute  in  Fairbanks. 

One  distinguishing  feature  of  the  Dali  City  earthquakes  (which 
they  were  named  after  an  abandoned  mining  camp  in  the  area)  was  that 
they  almost  invariable  showed  very  distinct  Pn  and  Pg  phases  at  nearly 
all  of  the  permanent  University  of  Alaska  stations  around  the 
epicentral  region.  In  addition,  an  intermediate  crustal  phase  was 
often  seen.  Depth  control  was  good  because  of  the  newly-installed 
temporary  stations,  and  it  was  recognized  that  this  presented  a  good 
opportunity  to  study  the  crustal  structure.  This  was  done,  and  the 
results  are  presented  in  the  following  section. 

STATUS  OF  RESEARCH 

The  Gold  King  Cluster 

In  the  search  for  earthquake  groupings  for  use  in  this  study,  it 
was  discovered  that  a  previously  unrecognized  cluster  could  be 
identified  near  a  landing  strip  at  Gold  King,  about  100  km  SSW  of 
Fairbanks.  These  earthquakes  nearly  all  occurred  within  a  volume  of 
about  20  km  in  diameter  (Figs.  1  and  2)  at  a  depth  centered  at  around 
100  km  (Fig.  3).  Although  the  first  event  which  could  be  identified 
as  belonging  to  the  group  was  recorded  in  1969,  most  of  the  analysis 
which  followed  used  only  data  dating  from  the  period  of  improved 
instrumental  coverage  from  1978  forward  (these  are  the  events  depicted 
in  Figs .  2  and  3 )  • 

As  cross-sections  A-A'  and  B-B'  show  (Figs.  1  and  3),  the  Gold 
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King  cluster  plainly  represents  the  northernmost  remnant  of  the 
subduction  zone  which  extends  southwestward  from  that  point  to  the 
western  edge  of  Cook  Inlet  and  out  along  the  Aleutian  Chain. 

A  composite  focal  mechanism  solution  obtained  from  37  of  the 
better-recorded  events  reveals  a  tensional  axis  dipping  steeply  to  the 
NNW  in  the  plane  of  the  subducting  plate.  This  indicates  that  these 
earthquakes  are  due  to  gravitational  sinking  within  the  plate.  (A 
more  detailed  presentation  will  appear  in  the  October,  1987,  issue  of 
the  Bulletin  of  the  Seismological  Society  of  America.) 

The  Pall  City  cluster 

Beginning  with  a  magnitude  (ML)  5-4  earthquake  on  February  14, 
1985,  a  three-month  series  of  strong  earthquakes  began  in  Central 
Alaska,  near  the  trans-Alaska  pipeline  about  30  km  north  of  the  Yukon 
River  crossing  of  the  pipeline  and  180  km  northwest  of  Fairbanks 
(Figs.  4  and  5).  This  series  culminated  in  a  magnitude  6.0  shock  on 
March  9.  which  was  the  largest  ever  recorded  this  far  north  in  the 
state.  In  all,  at  least  four  earthquakes  of  magnitude  5.0  or  better 
(Table  1),  and  a  total  of  38  events  of  magnitude  4.0  or  greater  were 
recorded  during  the  period. 

Because  this  is  an  unusual  area  for  significant  earthquakes  in 
Alaska,  the  series  aroused  the  immediate  interest  of  the  Geophysical 
Institute  seismology  staff.  Being  relative  accessible  (by  Alaskan 
standards),  it  was  soon  possible  to  establish  three  local  telemetered 
stations  around  the  site  (FISH,  RGV  and  BEH,  Fig.  5).  This  was  done 
despite  the  fact  that  it  was  mid-winter,  the  temperature  was  hovering 
at  around  -40,  and  there  were  40  knot  winds  at  times. 

The  development  of  the  cluster  with  time  can  be  seen  by  comparing 
the  two-to-three  week  interval  plots  in  Figure  6.  The  dispersion  in 
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Table  1 


largest  earthquakes  of  the  1985  north-central  Alaska  sequence  through 
March  16.  Coordinates  given  are  those  determined  at  the  Geophysical 
Institute.  The  magnitudes  shown  are  local  and  body-wave  magnitudes 
from  the  NOAA  Alaska  Tsunami  Warning  Center.  The  earthquake  at  111 6 
on  9  March  was  masked  on  Alaskan  records  by  preceding  events,  and  the 
jivsn  are  from  the  USGS  National  Sarthauake  Information  Service. 


DATE 

TIME(UT) 

MAGNITUDE 

DEPTH (KM) 

DAT ( N ) 

LONG ( W 

14  Feb  85 

0504 

5-^ 

13 

66.35 

149.82 

9  Mar  85 

1357 

4.8 

f 

0 

66 . 34 

149.83 

9  Mar  85 

1408 

6.0 

7 

66 . 35 

149.85 

9  Mar  85 

14 1 6 

5.* 

10 

66.2 

150.2 

9  Mar  85 

15^6 

1.5 

1 

66.33 

149.84 

9  Mar  85 

1621 

4.4 

13 

66 . 36 

149.87 

9  Mar  85 

1649 

4.3 

8 

66.3S 

149.87 

10  Mar  85 

0358 

4.4 

5 

66.31 

149.90 

10  Mar  85 

1210 

^.3 

16 

66.33 

149.91 

10  Mar  85 

1330 

5.6 

7 

66.30 

149.90 

10  Mar  85 

1916 

4.2 

8 

66.32 

149.9c 

10  Mar  85 

2305 

4.4 

2 

66.33 

149.9c 

16  Mar  85 

1333 

5.0 

17 

66.26 

149.33 
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Figure  6(A)  compared  with  the  relative  compactness  of  the  later  plots 
is  an  indication  of  the  refinement  in  location  capabilities  that  the 
installation  of  the  nearby  stations  permitted. 

Figure  7  demonstrates  that  the  Dali  City  earthquakes  were 
significantly  removed  from  activity  further  to  the  south,  although 
they  seem  to  lie  near  a  northern  extension  of  a  lineal  aftershock  zone 
that  followed  a  magnitude  6.5  earthquake  near  Rampart,  Alaska,  on 
October  9,  1968  (Gedney  et  al.,  1969;  Huang  and  Biswas,  1983). 

The  focal  mechanism  solution  obtained  is  best  interpreted  as 
left-lateral  faulting  on  a  NNE-SSW  trending  plane.  This  would  be  in 
accord  with  that  found  by  Huang  and  Biswas  ( 1 98 3 )  for  the  Rampart 
earthquake  to  the  south.  The  mechanism  for  the  March  9  event  also  has 
a  remarkable  parallel  in  an  earlier  event  in  the  same  area  exactly  1C 
years  earlier  on  March  9.  1975  (Estabrook  et  al.,  1986;  Fig.  8). 

With  the  nearby  stations  installed,  reasonably  accurate  focal 
depths  could  be  obtained.  To  take  advantage  of  at  least  two  very 
distinctive  crustal  phases  seen  on  the  seismograms  and  clear  Pn 
phases,  an  iterative  program  was  written  which  operated  on  the 
following  assumptions: 

(1)  Given  an  accurate  focal  depth  and  origin  time,  travel  times 
to  seismographic  stations  at  near-to-intermediate  (100  to  500  km) 
distances  are  determined  by  the  layer-velocity  structure  between 
source  and  receiver. 

(2)  Realizing  that  it  is  an  oversimplification,  if  one  assumes  a 
horizontally  plane-layered,  two-layer-over-half space  model,  the 
arrival  times  of  the  different  phases  at  individual  stations  at 
different  distances  should  be  a  function  of,  (a)  the  velocity  of  the 
basal  Pn  layer  at  the  Moho,  (b)  the  depth  of  this  layer,  (c)  the 
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velocity  of  layer  1  in  the  crust,  (d)  the  thickness  of  layer  1,  (e) 
the  velocity  of  layer  2  in  the  crust,  and  (f)  the  thickness  of  layer 

2. 

(3)  Factors  (a)  and  (c)  can  be  reasonably  well  determined  by 
observational  means,  so  that  the  observer  is  left  with  solving  for 
four  unknowns.  If  a  good  estimate  can  be  obtained  for  the  total 
crustal  thickness  (b) ,  solving  for  either  (d)  or  (f)  obviates  the  need 
for  determining  the  thickness  of  the  other  crustal  layer. 

(4)  The  consolidation  and  "interf ingering"  of  observations  of  Pn 
and  crustal  phases  at  a  number  of  different  locations  at  varying 
distances  around  the  epicentral  area  should  provide  an  indication  of 
which  crustal  model  best  matches  all  sites. 

To  find  the  best-fit  model,  an  exhaustive  program  was  written 
that  utilized  a  number  of  the  better-recorded  earthquakes  showing 
distinct  Pn  and  crustal  phases.  This  program  iterated  Moho  depths  and 
mantle  velocities  with  crustal  layer  thicknesses  and  crustal  layer 
velocities  in  every  plausible  combination.  It  then  matched  the  models 
with  the  observed  values  at  the  various  stations,  discarding  the 
"impossible”  models  and  retaining  those  which  might  have  produced  the 
observations.  The  result  was  an  enormous  volume  of  printout  (a  page 
of  which  is  reproduced  as  Figure  9)*  It  was  then  attempted  to  pick 
the  combination  which  best  matched  the  observed  observations  of  Pn  and 
the  crustal  phases  with  the  theoretical  calculations  at  as  many 
different  observing  stations  as  possible. 

In  actuality,  no  single  model  turned  out  to  be  "perfect."  This 
is  understandable  because  the  "ideal"  structural  model  does  not 
exists  absolutely  level,  plane  layering  in  the  Tanana-Yukon  Uplands 
cannot  be  expected,  although  the  study  area's  location  between  the 
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massifs  of  the  Brooks  and  Alaska  Ranges  is  probably  as  close  to  it  as 
can  be  expected  any  place  in  the  state. 

It  was  found  that  iterations  involving  increments  less  than  0.1 
km/sec  in  velocity  of  the  crustal  layers  or  1.0  km  in  thickness  were 
of  little  benefit.  Accordingly,  the  best-fit  parameters  for  the  area 
between  the  Dali  City  area  and  Fairbanks  were  found  to  be: 

Vp  (km/sec)  Depth  (km) 

6.0  0-24 

7.2  24-34 

7.91  +  0.12  below  34 

With  the  exception  of  the  depth  to  the  Moho,  these  values  agree 
well  with  those  obtained  by  Biswas  (unpublished)  which  have  been  used 
for  the  University's  routine  location  program  in  the  Alaskan  Interior. 
Biswas'  values,  obtained  largely  from  quarry  blasts  were: 

Vp  (km/sec)  Depth  (km) 

5.9  0-24.4 

7.4  24.4-40.2 

7.9  below  40.2 

It  should  be  noted  that  Biswas’  values  were  obtained  largely  from 
sites  further  to  the  south  and  closer  to  the  Alaska  Range.  It  might 
well  be  expected  that  greater  depths  to  the  Moho  would  be  measured 
there . 

The  Dali  City  sequence  also  provided  the  opportunity  to  make 
repeated  measurements  of  Pn  velocities  at  various  azimuths  radiating 
out  from  around  the  epicentral  area.  Referring  to  Figure  4,  these 
values  were: 
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Vp( km/sec) 

8.12  +  0.19 
8.08  +  0.29 
7.91  +  0.12 
7.65  +  0.14 
7.63  +  0.16 

In  this  table,  Azimuth  refers  to  the  direction  in  which  the  ray- 
leaves  the  epicenter  to  arrive  at  the  first  recording  station. 
Apparently,  the  variance  in  apparent  Pn  velocities  reflect  the 
downward  flexure  of  the  Moho  as  the  ray  approaches  the  Alaska  Range, 
slantingly  at  first,  and  then  more  at  right  angles. 

Supportive  of  this  concept  is  the  observation  that  the  leg  from 
RDS  to  WRH  is  7-91  km/sec  (which  has  here  been  taken  to  be 
representative  of  this  portion  of  the  Interior) ,  but  in  a 
straight-line  path  beyond  that,  as  the  Alaska  Range  is  approached  from 
WRH  to  GKC ,  the  apparent  velocity  drops  to  7.65  km/sec. 
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investigations  in  northern  Alaska"  was  based  largely  on  research  of 
the  Dali  City  earthquakes  described  above. 


Stations 

Azimuth 

GLM-HDA 

l4l 

FBA-HDA 

148 

RDS-WRH 

154 

WRH -GKC 

157 

NEA-MCK 

169 

8 


REFERENCES 


Estabrook,  C.,  J.  Davies,  L.  Gedney,  S.  Estes  and  J.  Dixon,  Seis 
micity  and  focal  mechanism  studies  of  the  1985  Dali  City 
earthquakes,  north-central  Alaska,  currently  undergoing 
review,  will  probably  appear  in  Bull.  Seism.  Soc.  Am.,  1986 

Gedney,  L.,  E.  Berg,  H.  Pulpan,  J.  Davies  and  W.  Feetham,  A 

field  report  on  the  Rampart,  Alaska  earthquake  of  October 
29,  1968,  Bull.  Seism.  Soc.  Am.,  59,  1421-1423,  1969. 

Huang,  P.  and  N.  Biswas,  Rampart  seismic  zone  of  central  Alaska, 
Bull.  Seism.  Soc.  Am.,  73,  813-829,  1983. 


9 


•  IMA 


50°  148°  146 


[  i  1  i  i  i 


FYU* 


TNN 


April  24,1974 


^  F®f#GLM 
/A  R£Stl  FAIRBANKS 
NEA  «CCB 
\  I  twrh^hda  / 


MCK 


AREA  ENLARGED  IN  FIG.  2 


Gold  King 
Cluster 


>DDM 


MCKINLEY 


-took- 


J  March  5, 

>  1973  %T0A 

•saw*scm 
•  PMR 


^ANCHORAGE/! 


i  i 


150° 


Prince 
^  William 

Sound 


Figure  1.  Three-letter  station  designators  show  the  locations  of  the  seismograp 
stations  utilized  in  this  study.  The  upper-hemisphere  focal  mechanisms  shown  all 
indicate  steeply-dipping  T-axes  (stippled  areas  indicate  congressional  quadrants) , 
Che  solution  labeled  Gold  King  cluster  is  a  composite  solution  obtained  during  this 
study.  Cross  section  A-A’  and  B-B'  are  shown  in  Figure  3  and  are  100  and  j0  km 
thick,  respectively. 


Figure  2.  Epicentral  plot  of  86  earthquakes  of  the  Gold  King  cluster  located  between 
1978  and  1984.  WRH,  GKC  and  MCK  are  local  seismographic  stations. 
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Figure  4.  (A)  Location  map  for  lower  part  of  figure.  Major  faults  and  lineaments  ar 
indicated.  (B)  The  Dali  City  area  enlarged  in  Figure  5  is  shown  as  the  stippled  area 
P5  and  P6  are  pump  stations  on  the  Trans-Alaska  pipeline. 
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Figure  5.  The  Dali  City  epicentral  area.  FISH,  RGV  and  BEH  are  three  temporary 
seismographic  stations  which  were  emplaced  shortly  after  the  earthquake  seouence  teg; 
The  location  of  the  main  shock  (star)  and  four  other  larger  shocks  are  indicated.  T’ 
tick  marks  on  the  streams  indicate  areas  of  broken  ice  and  stream  overflow.  The  cro 
hatched  area  indicates  a  hillside  along  which  several  springs  were  rejuvenated. 
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Figure  6.  Progression  of  epicentral  zone  development  from  February,  1985,  through 
April.  The  letters  indicate  focal  depths  in  10  km  increments,  i.e. ,  A=0-10  km 
B”10-20  km,  C=20-30  km,  etc.  ’ 
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Figure  9.  Sample  page  of  printout  from  iteration  program  which  matched  observed  Pn 
and  crustal  phases  with  various  possible  crustal  configurations.  Impossible 
combinations  were  not  printed 
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ABSTRACT 


Earthquake  data  collected  by  the  Geophysical  Institute 
of  the  University  of  Alaska - Fai rbanks  have  been  reprocessed 
to  study  the  seismicity  and  tectonics  of  northern  Alaska. 

A  seismic  velocity  structure  developed  in  the  Dali  City 
area  of  central  Alaska  suggests  granite  as  the  main 
component  of  the  upper  crust.  A  refraction  study  in  the 
eastern  Brooks  Range  constrains  the  Moho  dip. 
Microearthquakes  trend  east-west  roughly  parallel  to,  but 
south  of,  the  northern  foothills  of  the  Brooks  Range.  In 
northern  .'.laska  and  NW  Canada  the  Kobuk,  central  Kaltag, 
Porcupine,  Dali  Mt. -Rampart,  NW  Tintina  and  Richardson 
faults  are  seismically  active.  Eight  new  focal  mechanisms 
were  determined  for  the  northern  Alaska.  In  the  south  they 
show  strike-slip  faulting  consistent  with  the  historical 
movement  on  these  faults.  The  orientation  of  the  pressure 
axes  is  consistent  with  Pacific -North  American  plate 
convergence.  A  tectonic  model  developed  accounts  for  the 
pervasive  strike-slip  faulting. 
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Abstract.  Data  from  earthquakes  occurring  in 
northern  Alaska  collected  by  the  Geophysical  In¬ 
stitute  of  the  University  of  Alaska,  Fairbanks,  have 
been  reprocessed  using  a  seismic  velocity  model  de¬ 
veloped  for  the  Dali  City  and  Fort  Yukon  areas  of 
north  central  Alaska.  A  study  of  the  relocated  events 
shows  that  microearthquakes  occur  in  a  zone  roughly 
parallel  to,  but  south  of,  the  crest  of  the  Brooks 
Range.  The  events  also  show  that  the  Kobuk,  cen¬ 
tral  Kaltag,  Porcupine,  Dali;  Rampart,  northwest 
Tintina,  and  Eskimo  Lakes  faults  in  northern  Alaska 
and  northwest  Canada  are  seismically  active,  with 
the  Eskimo  Lakes  fault  being  most  active.  Eight 
new  focal  mechanisms  were  determined  for  northern 
Alaska  which  show  strike-slip  faulting  south  of  the 
Brooks  Range.  This  movement  is  consistent  with  the 
known  movement  on  the  associated  faults.  The  ori¬ 
entation  of  the  pressure  axes  derived  from  these  and 
other  fault  plane  solutions  is  consistent  with  Pacific- 
North  American  plate  convergence.  Large-scale 
shearing  appears  to  extend  northward  from  the 
Pacific-North  America  transform  boundary  into  arc¬ 
tic  Alaska  and  allows  the  interpretation  that  the 
Seward  Peninsula  extensional  zone  is  a  large-scale 
pull-apart  feature. 

Introduction 

Northern  Alaska  (Figure  1)  exhibits  a  higher  lev¬ 
el  of  seismicity  than  might  be  expected  considering 
the  distance  from  any  recognizable,  modern  plate 
boundary.  Several  moderate  to  large  earthquakes 
have  occurred  in  northern  Alaska:  a  magnitude  7.3 
PAS  (PAS  =  Ms  determined  by  the  California  Insti¬ 
tute  of  Technology)  event  in  the  Koyukuk  basin  in 
1958,  a  6.9  PAS  in  the  Chukchi  Sea  in  1928,  a  6.5  Ms 
(Ms,  Ml,  and  mb  a'  determined  by  the  National 
Earthquake  Informaaon  Service  (NEIS),  or  its  pre¬ 
decessors)  at  Rampart  in  1968,  a  5.2  Ml  (4.4  mb)  at 
Barter  Island  on  the  Arctic  coast  in  1968,  and  a  6.0 
Ms  at  Dali  City  in  the  Ray  Mountains  in  March  1985. 
To  help  understand  the. distribution  of  seismicity  and 
the  occurrence  of  moderate  to  large  earthquakes,  a 
crustal  velocity  model  has  been  determined.  Using 
this  model,  a  large  set  of  earthquakes  were  relocated, 
and  their  distribution  was  analyzed.  The  data  were 
also  used  to  derive  earthquake  fault  plane  solutions 
and  to  formulate  a  tectonic  model  to  explain  the  in¬ 
traplate  deformation  of  northern  Alaska. 

The  Brooks  Range  is  the  most  prominent  physio¬ 
graphic  feature  in  northern  Alaska  (Figure  1).  The 
range  trends  east-west  across  western  and  central 
Alaska  and  curves  toward  the  northeast  at  its  east- 
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ern  end.  Low-angle  north  verging  thrust  faults  with 
upward  of  500  km  of  lateral  shortening  characterize 
the  western  end  [Mayfield  et  al.,  1983),  but  these 
thrust  sheets  are  not  recognized  in  the  northeast 
Brooks  Range  but  may  be  present  farther  south 
[Oldow  et  al.,  1987).  Farther  to  the  east,  the  moun¬ 
tains  trend  south  along  the  west  side  of  the  Mac¬ 
kenzie  River,  eventually  becoming  the  Richardson 
Mountains,  which  are  separated  from  the  geological¬ 
ly  different  Brooks  Range  by  the  Porcupine  fault 
[Norris  and  Yorath,  1981]. 

The  area  south  of  the  Brooks  Range  is  dominated 
by  several  large  basins  and  a  bisecting  mountain 
range.  The  Yukon-Koyukuk  basin  contains  Creta 
ceous  sediments,  remnants  of  Mesozoic  island  arcs, 
and  is  intruded  by  Cretaceous  plutons  [Patton,  1973). 
The  basin  is  possibly  underlain  by  oceanic  crust  equi¬ 
valent  to  that  exposed  along  the  southern  edge  of  the 
Brooks  Range  and  along  the  western  flank  of  the 
Ruby  mountains.  To  the  east  of  the  Yukon-Koyukuk 
basin  is  a  large,  low-lying  area  of  Quaternary  sedi¬ 
ments  known  as  the  Yukon  Flats.  The  Yukon  Flats 
is  known  to  be  a  basin  of  considerable  depth,  but  the 
nature  of  its  basement  is  not  known.  Kirschner  et  al. 
[1985]  suggest  that  the  Yukon  Flats  developed  as  an 
extensional  graben  from  strike-slip  motion  along  the 
Tintina  fault,  though  most  accretion  models  imply 
that  it  may  be  underlain  by  oceanic  crust  (see,  for  in¬ 
stance,  Churkin  et  al.  [1982]). 

A  number  of  major  faults  have  been  mapped  in 
north  central  Alaska  (Figure  1),  but  the  relation¬ 
ships  between  them  are  often  obscure.  For  instance, 
the  Kaltag  fault,  which  trends  across  the  southern 
portion  of  the  study  area,  can  be  traced  northeast 
through  the  Yukon-Koyukuk  basin  and  Ruby  geanti¬ 
cline  before  disappearing  beneath  the  Yukon  Flats. 
Norris  and  Yorath  [1981]  believe  that  the  Kaltag 
fault  continues  northeast  beneath  the  Yukon  Flats  to 
emerge  as  the  Porcupine  lineament.  Another  possi¬ 
bility  is  that  it  may  wrap  around  the  southern  edge  of 
the  Y ukon  Flats  to  meld  with  the  Tintina  fault.  Sim¬ 
ilarly,  the  relationship  between  the  mapped  portions 
of  the  Porcupine  fault  and  the  Eskimo  Lakes  fault 
system  is  obscure,  as  is  the  connection  between  the 
Porcupine  fault  and  its  apparent  extension  as  the 
Porcupine  lineament. 

Crustal  Structure 

Previous  Studies 

There  are  no  published  crustal  models  that  apply 
specifically  to  northern  Alaska.  However,  several 
studies  have  been  made  in  adjacent  areas  (Table  1). 
Jin  and  Herrin  [1980]  studied  the  phase  and  group 
velocities  of  surface  waves  to  determine  an  average 
velocity  structure  for  south  central  Alaska.  Hanson 
et  al.  [1968]  derived  a  crustal  structure  model  from 
an  unreversed  refraction  line  across  the  Tanana 


12,026 


Estabrook  et  al.:  Seismotectonics  of  Northern  Alaska 


12,027 


Fig.  1.  Map  of  major  faults  in  Alaska.  BB,  Bruin  Bay;  BR,  Border  Ranges;  C,  Contact;  CH, 
Chiroskey;  CM,  Castle  Mountain-Lake  Clc.rke;  CS,  Chatham  Strait;  D,  Denali;  EL,  Eskimo 
Lakes;  F,  Fairweather;  FW,  Farewell;  HC,  Hines  Creek;  INF,  Innoko-Nixon  Fork;  K,  Kaltag; 
KT,  Kobuk  Trench;  P,  Porcupine;  Q,  Queen  Charlotte;  S,  Susitna;  T,  Tintina;  TC,  Totschunda. 
Also  shown  are  the  Yakutat  block  (YB)  and  other  key  locations  mentioned  in  the  text.  The 
shaded  areas  are  the  Yukon-Koyukuk  Province  and  the  Yukon  Flats.  Note  that  Dali  City  no 
longer  exists  but  is  used  for  the  general  location  of  the  Dali  City  earthquakes. 


Basin,  southwest  of  Fairbanks.  It  agrees  rather  well 
with  that  of  Jin  and  Herrin  [1980]  for  the  upper  15 
km.  Hanson  et  al.  ( 1968]  give  a  33-km  crustal  thick¬ 
ness  for  the  Fairbanks  area,  based  on  the  gravity 
models  of  Woollard  et  al.  [I960],  This  latter  thick¬ 
ness  is  considerably  shallower  than  the  43-km  Moho 
depth  of  Jin  and  Herrin  [1980],  which  is  an  average 
of  the  thick  crust  under  the  Alaska  Range  and  the 
thinner  surrounding  crust.  N.  Biswas  (University  of 
Alaska  Fairbanks,  unpublished  data,  1977)  deter¬ 
mined  a  crustal  velocity  structure  for  south  central 
Alaska  using  a  Monte  Carlo  inversion  of  earthquake 
and  quarry  blast  travel  times.  The  upper  crustal  ve¬ 
locity  of  5.9  km/s  for  the  top  24  km  is  slightly  faster 
than  that  of  Jin  and  Herrin  [1980]  and  Hanson  et  al. 
[1968].  The  crustal  thickness  of  Biswas's  model  is 
comparable  to  Jin  and  Herrin's  (40  and  43  km,  re¬ 
spectively),  but  thicker  than  Hanson  et  al.’s  [1968] 
for  the  Fairbanks  area  (33  km). 

Northern  Alaskan  Upper  Mantle  Velocity, 

Fort  Yukon  Study 

Refracted  phases  from  earthquakes  recorded  at 
stations  in  the  Barter  Island,  Fort  Yukon,  and  Fair¬ 
banks  networks  of  the  University  of  Alaska  (Figure 
2)  were  used  to  construct  a  reversed  seismic  refrac¬ 
tion  profile.  This  was  used  to  derive  an  upper  mantle 
or  Pn'velocity  under  the  Fort  Yukon  network  [Esta¬ 
brook,  1985]. 

A  well-located  earthquake  in  the  Fairbanks  area 
(781217;  Figures  2,  3a)  (all  events  are  labeled  as 
year,  month,  day,  e.g.,  781217  is  December  17,  1978) 


was  chosen  as  a  southern  source.  It  gives  a  Pn  veloci 
ty  of  7.94  ±0.18  km/s  (95%  confidence  limit  based  on 
work  by  Hanson  [1968]).  Earthquake  780621,  locat¬ 
ed  in  the  northeast  Brooks  Range  (Figures  2, 3b),  was 
chosen  as  a  northern  source.  This  event  yields  a  Pn 
velocity  of  7.88  ±0.22  km/s.  On  the  basis  of  this 
reversed  pair  of  values  (7.88  and  7.94  km/s)  across 
the  Fort  Yukon  network,  we  take  the  Pn  velocity  in 
northern  Alaska  to  be  7.9  km/s. 

Northern  Alaska  Crustal  Velocities.  Dali  City 
and  Fort  Yukon  Data 

The  close  proximity  of  seismic  stations  to  the  1985 
Dali  City  earthquakes,  events  7  and  8  in  Figure  8 
[Estabrook  et  al.,  1985],  provided  the  opportunity  to 
determine  a  local  P  and  S  wave  velocity  structure. 
To  do  this,  the  average  crusta1  thickness  between 
Dali  City  and  Fairbanks  was  fixed  at  35  km,  a  thick¬ 
ness  derived  by  L.  Gedney  (personal  communication, 
1986)  using  the  arrival  times  of  Pn  and  Py  waves. 
Wadati  diagrams  were  used  to  determine  the  Vp/Vs 
ratios  of  1.68 ±0.03  for  the  Dali  City  area.  This  is 
essentially  the  same  as  the  value  of  1.69  ±0.04 
determined  by  Gedney  and  Berg  [1969]  for  after 
shocks  of  the  1967  Fairbanks  earthquake.  By  sys 
tematically  varying  the  crustal  velocity  used  to  relo¬ 
cate  a  set  of  10  earthquakes,  it  was  possible  to  deter 
mine  the  velocity  for  which  hypocentral  and  travel 
time  errors  were  a  minimum.  These  errors  were 
large  for  velocities  less  than  6.0  km/s  and  asymptoti¬ 
cally  approached  a  small  lower  limit  for  faster  veloc¬ 
ities.  Thus  6.0  km/s  is  a  minimum  value  for  the 
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TABLE  1.  Velocity  Models  for  the  Crust  and  Upper  Mantle  of  Interior  and 

Northern  Alaska 


Number 

Model 

vP 

Depth 

Vp/V  s 

1 

Dali  City  area  (this  study) 

6.1 

0.0 

1.68 

7.9 

35.0 

1.68 

2 

above  model  with  mid-crustal  layer 

6.1 

0.0 

1.68 

6.6 

20.0 

1.68 

7.9 

35.0 

1.68 

Fairbanks  Area 

3 

upper  crust  from  Hanson  [1988]; 

4.5 

0.0 

1.68 

lower  crust  and  Moho  from 

5.4 

2.0 

1.68 

refraction  in  Fort  Y ukon  area 

5.8 

4.5 

1.68 

(this  study) 

6.1 

11.0 

1.68 

8.0 

33.0 

1.68 

Brooks  Range 

5.4 

0.0 

L.68 

5.8 

4.5 

1.68 

6.1 

11.0 

1.68 

8.0 

40.0 

1.68 

4 

standard  interior  Alaska  model  (N. 

5.9 

0.0 

1.78 

N.  Biswas,  unpublished  data. 

7.4 

24.4 

1.78 

1977) 

7.9 

40.2 

1.78 

5 

above  model  with  Vp/Vs  ratio  from 

5.9 

0.0 

1.68 

this  study 

7.4 

24.4 

1.68 

7.9 

40.2 

1.68 

Vp,  velocity  of  p  wave  in  km/s.  Depth,  depth  of  the  top  of  the  layer  in  km.  Vp / 
Vs>  ratio  of  the  velocity  of  the  p  wave  to  that  of  the  s  wave;  used  in  the  hypo- 
center  location  program  to  convert  s  wave  arrival  times  to  equivalent  p  wave 
times. 


average  velocity  of  the  upper  10  km  of  the  crust  in 
this  area. 

Crustal  velocities  were  also  determined  using 
events  located  within  and  recorded  by  the  Fort 
Yukon  network  (Figure  2).  Figure  3c  shows  a  good 
fit  to  a  straight-line  relationship  between  arrival 
time  and  distance,  even  though  the  stations  surround 
the  epicenter.  This  is  interpreted  as  representing  a 
velocity  of  6.04±0.12  km/s  beneath  the  network.  On 
the  basis  of  the  Dali  City  and  Fort  Yukon  data,  a  P 
wave  velocity  of  6.1  km/s  is  taken  for  the  upper  crust 
of  northern  Alaska. 

Velocity  Model  Tests 

The  combined  results  of  the  Fort  Yukon  and  Dali 
City  studies  yield  a  P  wave  velocity  model  consisting 
of  a  6.1  km/s  crust,  35  km  thick,  overlying  a  7.9  km/s 
upper  mantle.  Location  errors  and  travel  time  resid¬ 
uals  for  this  model  are  compared  in  Table  2  to  the 
models  of  Hanson  et  al.  [1968]  and  N.  Biswas 
(University  of  Alaska  Fairbanks,  unpublished  data, 
1977). 

When  the  entire  northern  Alaska  data  set  is  con¬ 
sidered  and  only  a  sparse  regional  network  is  used, 
the  similarities  in  the  estimates  of  error  from  one 
model  to  the  other  indicate  that  there  is  no  obvious 
preferred  model.  However,  when  local  stations  are 
involved,  as  shown  for  Dali  city  in  Table  2,  the  errors 
increase  for  models  with  lower  velocities  in  the  crust 


(models  3,  4,  and  5).  Therefore  velocity  model  1,  Ta 
ble  1,  has  been  used  for  this  study  since  it  has  com¬ 
parable  errors  to  other  models  on  a  regional  scale  (at 
stations  farther  than  the  Pn  crossover  distance)  and 
has  the  advantage  of  simplicity. 

Seismicity 

The  present  study  incorporates  all  of  the  data  col¬ 
lected  by  networks  operated  by  the  Geophysical  In¬ 
stitute  of  the  University  of  Alaska  for  the  years  1967- 
1984. 

The  set  of  earthquake  phase  data  for  these  years 
was  compiled  into  a  single  file  (5427  events),  and  the 
earthquakes  were  relocated  using  HYPOELLIPSE 
[Lahr,  1982]  with  a  single-layer  crustal  model  (model 
1,  Table  1). 

Using  the  selection  criteria  that  the  horizontal  lo¬ 
cation  errors  (ERH)  be  less  than  10  km,  the  depth  er 
rors  (ERZ)  less  than  20  km,  and  the  root  mean  square 
of  the  travel  time  residuals  less  than  1  s,  a  subset  was 
created  which  sharpened  many  clusters  and  trends  of 
earthquakes  (Figure  4).  Note  that  the  detection 
threshold  is  not  homogeneous  in  time  nor  space  for 
the  earthquakes  whose  epicenters  are  shown  in  Fig 
ure  4.  A  more  homogeneous  representation  of  the 
seismicity  of  northern  Alaska  is  given  in  Figure  5 
which  shows  an  epicenter  map  of  all  earthquakes 
greater  than  magnitude  3  during  the  years  1976 
1979.  During  this  period  the  Seward  Peninsula,  Fort 
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Fig.  2.  Location  map  and  station  codes  for  seismo¬ 
graph  stations  in  interior  and  northern  Alaska.  The 
Seward,  Barter  Island  (BI),  and  Fort  Yukon  (FY) 
nets  were  operating  together  from  1976  to  1979.  The 
refraction  study  used  earthquakes  (circles)  designat¬ 
ed  by  year,  month,  day,  and  all  the  stations  east  of 
155°W  with  the  exception  of  BM3. 


Yukon,  and  Barter  Island  networks  were  operating 
simultaneously,  and  with  the  exception  of  a  small  re 
gion  near  Point  Barrow,  the  detection  threshold  of 
the  combined  network  is  estimated  to  be  about  mag¬ 
nitude  3.0  [Estabrook,  1985].  This  estimate  is  based 
largely  on  magnitude-distance  plots  for  each  station. 
Figure  5  shows  a  rather  uniform  distribution  of 
events  along  the  length  of  the  Brooks  Range  south  of 
the  continental  divide.  The  lack  of  magnitude  3  and 
larger  events  north  of  the  Brooks  Range  on  the  Arctic 
Slope  during  the  1976-1979  time  window  may  not  be 
representative  of  the  long-term  level  of  seismicity 
since  the  time  period  is  so  short. 

However,  the  same  absence  of  events  on  the  Arc¬ 
tic  Slope  is  seen  in  Figure  6,  which  shows  epicenters 
for  earthquakes  of  magnitude  4  and  larger.  This 
map,  based  on  data  from  the  National  Earthquake 
Information  Service  [1967-1985],  is  probably  com¬ 
plete  since  1970  for  magnitude  4  and  since  1960  for 
magnitude  5  and  larger  events.  A  clear,  more  or  less 
north-south  trend  of  these  events  near  the  Alaska- 
Canada  border  can  be  seen  in  Figure  6.  An  appar¬ 
ently  separate  distribution  of  larger  events  shows  a 
distinct  west-northwest  trend  from  the  vicinity  of 
Fairbanks,  as  noted  by  Davies  [1983], 

In  addition  to  the  observation  that  most  of  the 
events  lie  to  the  south  of  the  continental  divide,  there 
are  several  clusters  and  trends  that  can  be  seen  on 
the  maps  of  earthquakes  north  of  65°N.  These  fea¬ 
tures  are  numbered  in  Figure  4  and  described  in  Ta¬ 
ble  3. 

Focal  Mechanisms 

Eight  new  fault  plane  solutions  were  determined 
for  earthquakes  in  northern  Alaska  (Figure  7). 
These  solutions  are  composed  of  P  wave  first  motions 
plotted  on  equal-area  lower  hemisphere  projections 
ef  the  focal  sphere.  The  first  motions  are  predomi¬ 
nantly  from  Alaskan  regional  stations  and  from  sta¬ 
tions  of  the  Canadian  Seismograph  Network.  For  the 
larger  events  (magnitude  greater  than  4.5),  stations 
from  the  World-Wide  Standard  Seismograph  Net¬ 
work  were  also  used.  The  Alaskan  stations  were 


checked  for  polarity  ambiguities  using  nuclear  explo¬ 
sions  and  teleseisms  as  calibration  events  [Esta- 
brook,  1985]. 

The  depths  of  the  earthquakes  whose  focal  mecha 
nisms  were  determined  in  this  study  were  fixed  at  10 
km  because  of  poor  depth  control.  This  depth  was 
chosen  because  results  from  studies  of  the  Dali  City 
earthquakes  [Estabrook  et  al.,  1985]  and  of  earth 
quakes  in  the  Fairbanks  area  [Gedney  et  al.,  1980] 
show  that  most  interior  and  northern  Alaska  earth 
quakes  occur  within  the  upper  20  km  of  the  crust. 
Model  1  (Table  1)  was  used  to  calculate  takeoff 
angles  and  to  locate  the  earthquakes. 

Stations  that  are  farther  away  than  the  crust- 
mantle  crossover  distance  (about  150  km)  are  plotted 
at  a  common  takeoff  angle  regardless  of  their  actual 
epicentral  distance.  This  condition  holds  until 
roughly  1500  km  (about  15°),  at  which  distance  deep- 


100  200  300  *00  500 

DISTANCE  (km) 

Fig.  3.  Reduced  travel  time  versus  distance  plot  for 
earthquakes (a) 781217, (b)  780621,  and  (c) 770825. 
Larger  circles  denote  station  arrivals  used  in  the 
velocity  determination;  7.6  km/s  is  the  reduction  ve¬ 
locity. 
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TABLE  2.  Velocity  Model 


Entire  Northern 
Alaska  Data  Set 


Only  Regional  Stations 
at  Dali  City 


Local  and  Regional 
Stations  at  Dali  City 


Number  Number  Number 

of  of  of 


Model 

Events 

Mean 

S.D. 

Events 

Mean 

S.D. 

Events 

Mean 

S.D. 

1  Depth 

5427 

16.16 

31.57 

20 

7.79 

7.31 

47 

13.49 

3.66 

rms 

5427 

0.58 

1.96 

20 

0.60 

0.31 

47 

0.20 

0.19 

ERH 

5427 

30.37 

37.35 

20 

9.52 

21.02 

47 

3.33 

1.78 

ERZ 

5427 

48.66 

44.06 

20 

40.64 

43.12 

47 

4.28 

2.49 

2  Depth 

5423 

15.80 

31.13 

20 

10.01 

14.37 

47 

11.75 

4.29 

rms 

5423 

0.59 

1.42 

20 

0.64 

0.34 

47 

0.20 

0.20 

ERH 

5423 

31.79 

37.71 

20 

8.67 

21.12 

47 

3.32 

1.79 

ERZ 

5423 

48.78 

44.01 

20 

45.61 

44.33 

47 

4.42 

2.38 

3  Depth 

5430 

14.04 

20.44 

19 

8.24 

8.12 

47 

7.74 

3.97 

rms 

5430 

0.60 

1.29 

19 

0.61 

0.32 

47 

0.23 

0.18 

ERH 

5430 

34.11 

38.57 

19 

10.95 

21.81 

47 

6.81 

7.82 

ERZ 

5430 

45.09 

42.31 

19 

39.12 

42.52 

47 

30.00 

35.59 

4  Depth 

5426 

16.95 

25.54 

20 

27.28 

17.32 

47 

5.14 

4.33 

rms 

5426 

0.58 

1.28 

20 

0.52 

0.28 

47 

0.25 

0.14 

ERH 

5426 

33.19 

38.29 

20 

8.08 

21.27 

47 

2.68 

1.34 

ERZ 

5426 

47.45 

42.85 

20 

35.02 

38.61 

47 

32.41 

38.71 

5  Depth 

5427 

15.83 

26.39 

20 

12.62 

15.08 

47 

9.31 

4.59 

rms 

5427 

0.58 

1.28 

20 

0.64 

0.34 

47 

0.21 

0.18 

ERH 

5427 

34.45 

38.55 

20 

10.84 

23.39 

47 

3.31 

2.29 

ERZ 

5427 

45.23 

42.30 

20 

23.22 

33.08 

47 

6.63 

6.56 

Statistical  tests  of  the  velocity  models  in  Table  1.  In  these  tests  a  given  set  of  earthquakes  are  relocated 
using  each  of  the  five  models,  and  the  mean  depths  and  standard  deviations  are  compared.  Also  compared 
are  the  means  and  standard  deviations  of  three  measures  of  the  quality  of  the  hypocenter  locations:  (1)  the 
root-mean-square  difference  between  the  model  and  observed  arrival  times  (rms);  (2)  and  (3)  the  horizontal 
and  vertical  projections,  ERH  and  ERZ  of  the  error  ellipse  calculated  by  HYPOELLIPSE  [Lahr,  1982], 
Three  different  data  sets  are  used:  (1)  the  entire  northern  Alaska  data  set.  some  5430  events:  (2)  a  set  of  20 
events  in  the  Dali  City  area  located  using  only  regional  stations  (no  local  stations  of  Dali  City);  and  (3)  a  set 
of  47  events  located  using  both  local  and  regional  stations.  Significant  differences  appear  for  (1)  the  mean 
depth  calculated  using  model  4  for  the  regional-only  data  set;  (2)  and  (3)  the  mean  ERH  and  standard 
deviations  for  models  3  and  4,  respectively,  for  the  local  and  regional  data  set.  S.D.,  standard  deviation. 


er  layers  of  the  Earth  begin  to  be  sampled.  Takeoff 
angles  for  stations  farther  than  1500  km  are  taken 
from  the  Herrin  tables  [Herrin  et  al.,  1968). 

Earthquake  fault  plane  solutions  have  been  com¬ 
piled  for  three  different  regions  (Figure  8  and  Table 
4):  (1)  east  central  Interior  and  the  Yukon  Territory 
(events  1-9);  (2)  western  Alaska,  including  parts  of 
the  Yukon-Koyukuk  basin,  Seward  Peninsula,  and 
the  western  Brooks  Range  (events  10-16);  and  (3) 
northeast  Alaska  and  the  Canadian  margin  (events 
17-23). 

Region  1:  East  Central  Interior  Alaska 
and  Yukon  Territories 

Mechanisms  for  region  1  are  all  very  similar:  all 
are  strike-slip  solutions  with  northwesterly  trending 
pressure  axes.  Event  1  (this  study)  is  a  Ms  5.0  earth¬ 
quake  on  the  Tintina  fault  on  the  south  side  of  the 
Yukon  Flats.  Here  the  east-west  trending  nodal 
plane  is  chosen  as  the  fault  plane,  indicating  right- 
lateral  strike-slip  motion.  This  is  consistent  with  the 
historical  movement  on  this  major  fault  [Gabrielse, 
1985].  Event  2  was  a  strike-slip  mechanism  (Cook  et 


al.,  1984;  Coley,  1983],  The  northeast  trending  nodal 
plane  is  chosen  as  the  fault  plane  because  it  parallels 
mapped  faults  as  shown  by  Beikman  [1980], 

Events  3  [Gedney,  1970]  and  4  [Jordan  et  al., 
1968]  reflect  right-lateral  movement  on  northwest 
trending  fault  planes  which  parallel  the  trend  of  the 
Fairbanks  seismic  zone  [Gedney  and  Berg,  1969; 
Gedney  et  al.,  1980].  Event  5,  the  Ms  6.5  Rampart 
earthquake  [Gedney  et  al.,  1969;  Huang,  1979; 
Huang  and  Biswas,  1983;  Coley,  1983],  has  a  strike- 
slip  mechanism  similar  to  the  Dali  City  (events  7  and 
8)  and  Fairbanks  area  events.  Huang  [  1979],  Huang 
and  Biswas  [1983],  and  Coley  [1983]  suggested  that 
the  faulting  of  the  Rampart  event  was  left-lateral 
strike-slip  on  a  north-south  trending  plane  because  it 
paralleled  the  aftershock  zone  (Figure  7). 

Events  6,  7,  and  8  (this  study)  have  left-lateral 
strike-slip  mechanisms  on  north  to  northeast  trend¬ 
ing  fault  planes.  The  north-south  trending  nodal 
planes  were  chosen  on  the  basis  of  a  detailed  study  ol 
the  Dali  City  earthquake  sequence  which  reveals 
that  the  largest  events  (magnitude  greater  than  4.5) 
occurred  near  the  middle  of  a  north-northeast  trend 
ing  elongate  aftershock  zone  (Estabrook  et  al..  19851 


Estabrooketal.:  Seismotectonics  of  Northern  Alaska 


12,031 


Fig.  4.  Earthquake  epicenters  for  the  period  1967-1984  with  horizontal  location  errors  (ERH) 
<  10  km,  depth  errors  (ERZ)  <20  km,  root-mean-square  travel  time  residuals  (RMS)  <  1.00  s 
and  number  of  stations  >4.  Dashed  line  represents  the  continental  divide.  Numbers  refer  to 
clusters  of  events  described  in  Table  3. 


and  close  to  the  active  Dali  Mountain  fault  (which 
trends  north-south)  [Brogan  et  al.,  1975]  (G.  Plafker, 
unpublished  data,  1986).  Event  6,  located  to  the 
south  of  the  other  two  in  this  group,  indicates  a  pos¬ 
sible  southward  continuation  of  the  Dali  Mountain 
fault.  A  minor  additional  aftershock  trend  coincides 
with  an  east-northeast/south-southwest  trending 
topographic  high,  which  could  be  the  surface  expres¬ 
sion  of  an  unmapped  fault  subparallel  to  the  Kaltag 
fault  and  Porcupine  lineament. 

Event  9  (this  study),  a  Mb  4.6  earthquake  on  the 
Kobuk  Trench,  has  also  been  studied  by  Gedney  and 
Marshall  [1981],  It  is  grouped  with  other  events 
from  interior  Alaska  because  they  all  have  similarly 


oriented  strike-slip  solutions.  The  solution  from  this 
study  has  nodal  planes  with  slightly  different  orien 
tations  than  those  of  Gedney  and  Marshall,  as  a  re 
suit  of  incorporating  data  from  the  International 
Seismological  Centre  (ISC)  [1983]  bulletin.  Critical 
to  the  determination  of  the  strike  of  the  east-west 
nodal  plane  are  the  regional  Alaskan  stations  FYU 
and  BM3  (Figure  2).  These  two  stations  have  almost 
identical  waveforms,  yet  opposite  first  motions.  The 
station  polarities  were  thoroughly  checked  and  could 
not  explain  the  discrepancy.  The  only  other  possibil¬ 
ity  is  that  of  a  polarity  changing  reflection  along  the 
ray  path.  The  east-west  trending  plane  is  chosen  as 
the  fault  plane  because  it  is  parallel  to  the  strike  of 


Fig.  5.  Earthquakes  with  magnitudes  greater  than  3  for  the  years  1976-1979  when  Seward 
Peninsula,  Barter  Island,  and  Fort  Yukon  nets  were  operational.  When  these  networks  were 
operating,  the  magnitude  threshold  for  much  of  northernmost  Alaska  was  smaller  than  mag 
nitude  3.  Dashed  line  represents  the  continental  divide. 
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Fig.  6.  Earthquakes  with  magnitudes  greater  than  4  for  all  years.  Source  NEIS.  Faults  are 
from  King  [1969]. 


the  presently  active  Kobuk  Trench  [Brogan  et  al., 
1975], 

Region  2;  Western  Alaska 

Earthquake  focal  mechanisms  from  western 
Alaska  (events  10-16)  can  be  divided  into  two  groups; 
those  north  of  66°N  are  strike-slip,  while  those  south 
of  that  latitude  give  normal  faulting  solutions. 

Solutions  by  several  workers  [Ritsema.  1962; 
Fara,  1964]  for  event  10,  the  1958  Huslia  earthquake 
(Ms  7.3),  give  a  normal  fault  mechanism  on  a  north¬ 
east  trending  plane.  Davis  [1960]  mapped  ground 
breakage  with  a  similar  trend.  A  normal  fault  mech¬ 
anism  is  atypical  for  this  area.  Event  11  (this  study) 
is  a  composite  fault  plane  solution  determined  for 
four  events  with  magnitudes  greater  than  4.0  that 
occurred  on  the  Purcell  Mountain  trend.  If  this  trend 
is  the  fault  plane,  then  event  11  represents  a  left- 
lateral  strike-slip  fault.  It  is  not  clear  whether  the 
Purcell  Mountain  trend  is  related  to  event  10,  the 
Huslia  earthquake. 

Event  12,  a  mb  5.2  earthquake  which  occurred  in 
1981  north  of  Kotzebue,  has  been  studied  by  D.  B. 
Cook  (personal  communication,  1985),  Coley  [1983], 
and  Biswas  etal.  [1983].  Cook  and  Coley  selected  the 
east-west  trending  nodal  plane  as  the  probable  fault 
plane  because  it  is  parallel  to  the  trend  of  aftershocks 
which  they  relocated  using  the  master  event  tech¬ 
nique.  This  solution  is  used  here  because  the  solu¬ 
tion  to  Biswas  et  al.  [1983]  is  a  composite  with  many 
inconsistent  readings.  Event  13  [Coley,  1983]  is  a 
strike-slip  solution  similar  to  event  12,  possibly  in¬ 
dicating  right-lateral  motion  on  the  Kobuk  Trench. 
These  mechanisms  have  nodal  plane  orientations 
similar  to  other  strike-slip  events  in  central  interior 
Alaska. 

Earthquakes  14,  15  (D.  B.  Cook,  personal  com¬ 
munication,  1985),  and  16  [Liu  and  Kanamori,  1980; 
Sykes  and  Sbar,  1974]  in  the  southern  Seward  Penin¬ 
sula  and  Norton  Sound  area  show  mainly  normal 
faulting  on  east-west  oriented  planes.  This  agrees 
with  the  general  trend  of  mapped  Neogene  normal 
faults  [Hudson  and  Plafker,  1978],  Event  14  (D.  B. 
Cook,  personal  communication,  1985)  is  a  normal 


fault  with  a  strike-slip  component,  with  nodal  plan' 
perpendicular  to  events  12  and  13. 

Region  3:  Northeast  Alaska  and  Beaufort  Sea 

A  focal  mechanism  was  determined  for  event 
(this  study),  an  earthquake  along  the  continent 
margin  of  the  Beaufort  Sea.  This  complements  tl 
solutions  determined  by  Fujita  et  al.  [1983]  for  tl 
Barter  Island  Mb  4.7  earthquake  (event  18)  and  tv 
from  the  Canadian  margin  (event  20  and  21)  dete 
mined  by  Hasegawa  et  al.  [1979]  and  Hasegav 
[1977],  respectively.  A  somewhat  different  solutii 
was  obtained  by  Biswas  et  al.  [1986]  for  event  1 
however,  the  solution  of  Fujita  et  al.  [1983]  was  us. 
because  of  the  verification  of  the  solution  by  surfa 
wave  radiation  patterns.  All  of  the  solutions  fro 
along  the  margin  (events  18,  19,  20,  and  21)  she 
significant  strike-slip  components  but  are  incons. 
tent  in  nodal  plane  orientation.  Event  20,  an  mb  £ 
in  the  Mackenzie  Delta  area,  has  a  P  axis  orient 
perpendicular  to  the  spreading  direction  of  the  arc 
spreading  center.  Event  19  (this  study)  is  similar 
event  20  [Hasegawa  et  al.,  1979],  possibly  indicatii 
a  similar  tectonic  process  in  the  region  of  the  Mt 
kenzie  River  delta.  However,  the  P  axes  of  events 
(D.  B.  Cook,  personal  communication,  1985)  and 
[Hasegawa  et  al.,  1979]  are  rotated  90°  compar 
with  those  of  events  19  and  20. 

The  solution  for  event  17  (this  study)  is  either 
normal  faulting  mechanism  or  a  thrust  mechanisi 
but  in  either  case  the  tension  axis  trends  nort 
south.  The  normal  fault  solution  is  preferred  becau 
the  east-west  trending  nodal  plane  parallels  t 
structural  grain  of  the  northeast  Brooks  Ranf 
These  regionally  inconsistent  solutions  cannot 
reconciled  unless  one  invokes  slip  along  preexisti. 
fractures  or  local  stress  fields  acting  in  place  of  t 
regional  one  for  some  of  the  events  discussed. 

Two  solutions  have  been  determined  for  eart 
quakes  in  the  Richardson  Mountains  of  the  Yuk 
Territory.  Event  22  is  a  very  well  constrained  stril 
slip  solution  (LeBlanc  and  Wetmiller,  1974]  whi 
falls  on  the  north-south  trending  Eskimo  Lakes  fai 
system.  If  the  fault  plane  is  north-south,  then  t 


Estabrook  et  al.:  Seismotectonics  of  Northern  Alaska 


12,033 


TABLE  3.  Earthquake  Clusters 


No. 

Name 

Location 

Description 

1 

Ear  Mountain 
cluster 

65.7°N 

166. 0°W 
Figure  4 

Seward  Peninsula:  The  largest.  Ml  2.7,  and  the  other  27  events, 
were  recorded  throughout  the  lifetime  of  the  Seward  Penin¬ 
sula  network  ( 1977- 1983). 

2 

Baird  Mt.  cluster 

67.8°N 

161. 5°W 
Figure  4 

North  of  Kotzebue  main  shock  (mb  5.2  1982)  and  aftershocks. 
This  cluster  contains  12  events  with  mb  greater  than  4.0. 

3 

Purcell  Mountains 
trend 

66.5°N 

157. 5°W 
Figure  4 

North-south  trend  of  86  events,  five  have  mb  greater  than  4.0. 
Activity  pronounced  between  1971  and  1973,  though  earth 
quakes  continue  to  the  present. 

4 

Indian  Mountain 
cluster 

66.0°N 
155.5°W 
Figure  4 

Active  since  the  beginning  of  seismic  recording  by  the  Geophysi 
cal  Institute.  Most  activity  and  largest  event.  Ml  3.5,  occur 
red  in  1976.  Many  events  with  Ml  between  2  and  3. 

5 

Kobuk  trench 

67.0°N 

154. 7"W 
Figure  4 

A  1980  sequence  of  magnitude  4-5  earthquakes  [Gedney  and 
Marshall,  1981],  Numerous  aftershocks  were  recorded  at 
IMA,  the  closest  station. 

6 

Rampart  trend 

65.5°N 
150.0°W 
Figure  4 

1968  Rampart  sequence,  a  1968  magnitude  6.5  event  and  its 
north-south  trending  aftershocks,  120  km  northwest  of  Fair¬ 
banks.  A  resurgence  of  activity  occurred  in  1972  [see  Huang, 
1979;  Huang  and  Biswas,  1983], 

7 

Dali  City  after 
shocks 

66.25°N 
150. 0°W 
Figure  6 

These  earthquakes,  events  7  and  8  in  Figure  8,  were  recorded 
between  February  and  April  1985  (continuing  to  the  present) 
and  thus  do  not  appear  in  Figure  4;  the  largest  was  a  magni 
tude  6.0  (Ms)  (Estabrook  etal.,  1985], 

8 

Chandalar  River 
cluster 

66.9°N 
146.3°W 
Figure  4 

Ml  3.5  event  occurred  in  a  sequence  of  14  events  during 
January-March  1979.  Total  of  19  events.  It  is  noteworthy 
that  this  cluster  and  the  Porcupine  River  cluster  were  active 
at  nearly  the  same  time. 

9 

Tintina  fault 

65.7°N 
145. 2°W 
Figure  4 

A  line  of  events  that  coincides  with  the  Tintina  fault.  Much  of 
this  activity  is  in  the  Tintina  Trench,  northwest  of  central 
Alaska. 

10 

Porcupine  River 
cluster 

66.9°N 
144.0°W 
Figure  4 

Mi,  3.0  is  longest  event  in  a  cluster  of  11  events,  mostly  during 
November  1978  to  January  1979,  located  on  the  mapped  trace 
of  the  Porcupine  lineament. 

11 

Sagavanirktok 
River  cluster 

68.5°N 
149.0^ 
Figure  4 

Main  shock  and  aftershock  sequence  of  a  mb  4.3  earthquake. 
May  1980. 

12 

Barter  Island 
swarm 

71.3°N 
143. 5°W 
Figure  8 

Barter  Island  swarm,  January  1968,  largest  event  Ml  5.2,  occur 
red  prior  to  the  installation  of  northern  Alaska  stations.  The 
Canadian  Seismograph  Network  recorded  five  events  with 
magnitudes  greater  than  4.0  in  this  swarm. 

13 

Eskimo  Lakes 
fault  (also 
known  as  the 
Richardson 
fault) 

66.0°N 
135. 0°W 
Figures  1 
and  8 

Historical  earthquake  file  of  NEIS  gives  five  events  with  magni 
tudes  greater  than  6.0  and  109  events  with  magnitudes  great 
er  than  4.0.  Earthquakes  associated  with  the  Richardson 
Mountains  and  the  Rocky  Mt.  Front.  A  temporary  local  net 
work  and  determined  that  the  seismicity  was  associated  with 
a  north-south  trending  right-lateral  strike-slip  fault  [LeBlanc 
and  Wetmiller,  1974], 

focal  mechanism  is  right-lateral  strike  slip.  Event  23 
(Sykes  and  Sbar,  1974)  is  an  ambiguous  thrust  mech¬ 
anism,  yet  the  fault  planes  are  parallel  to  mapped 
thrust  faults  in  the  Canadian  Cordillera  (King, 
1969).  These  two  events  have  P  axes  oriented  in  a 
northerly  direction  consistent  with  a  common  tecton 
ic  stress  system. 


Discussion 

The  seismicity  distribution  and  the  focal  mecha 
nisms  determined  in  this  study  and  by  others,  togeth¬ 
er  indicate  that  pervasive  intraplate  deformation  is 
occurring  in  Alaska  and  northwest  Canada,  well  to 
the  north  of  the  Alaskan  subduction  zone.  This  study 
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Fig.  7.  Earthquake  fault  plane  solutions  determined  for  this  study.  Triangles  are  compres- 
sional  (solid)  and  dilatational  (open)  first  motions.  Small  triangles  are  station  readings  from 
the  ISC  catalog.  The  letters  "P”  and  "T”  are  the  pressure  and  tensional  axes,  respectively. 


has  also  shown  or  confirmed  that  many  of  the  major 
crustal  faults  in  northern  Alaska  and  northwest 
Canada  are  seismically  active.  Among  these  are  the 
Tintina,  Kaltag,  Kobuk,  Eskimo  Lakes  faults,  and 
possibly  the  Porcupine  lineament  (Figure  1).  Fur¬ 
thermore,  fault  plane  solutions  for  earthquakes  on 
the  Tintina,  Kobuk,  and  Eskimo  Lakes  faults  reveal 
that  the  present  faulting  is  in  the  same  sense  as  their 
historical  movement.  However,  as  the  seismicity 
maps  reveal,  most  of  the  earthquakes  smaller  than 
magnitude  5  are  not  aligned  along  these  major  faults 
and  lineaments. 

A  number  of  workers  [e.g.,  Jones  et  al.,  1981: 
Churkin  et  al.,  1982;  Moore  et  al.,  19851  have  sug 


gested  that  central  and  northern  Alaska  are  cor  | 
posed  of  tectonostratigraphic  terranes,  as  is  southe' 
Alaska.  By  definition,  these  terranes  are  fai 
bounded,  and  the  seismicity  associated  with  some 
these  deep  crustal  faults  may  indicate  that  these  tf 
ranes  are  still  responding  to  externally  applr 
stresses.  The  seismicity  within  the  terranes  is  mo  ! 
enigmatic,  possibly  occurring  on  reactivated  min 
faults  and  other  zones  of  weakness.  The  largt 
earthquakes  in  northern  Alaska  do  not  always  occ 
at  terrane  boundaries  or  along  a  single  major  fau 
Of  the  earthquakes  larger  than  magnitude  6,  eari 
quakes  along  the  Eskimo  Lakes  fault  (LeBlanc  a 
Wetmiller,  1974]  and  events  on  the  Tintina  fault  f 
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Fig.  8.  Earthquake  fault  plane  solutions,  magnitude  4  and  larger  earthquakes  (NEIS),  and 
geologic  faults  [King,  19691.  Fault  plane  solutions  are  equal-area  lower  hemisphere  projec¬ 
tions  showing  quadrants  of  compressional  (shaded)  and  tensional  first  motion.  Mechanisms 
are  oriented  with  north-south  parallel  to  the  edge  of  the  paper,  not  parallel  to  map  north. 
Mechanism  parameters  are  found  in  Table  4. 


associated  with  major  tectonic  boundaries.  Other 
large  earthquakes  (Huslia,  Fairbanks,  and  Dali  City) 
were  not  on  known  terrane-bounding  faults,  which 
supports  the  suggestion  that  there  must  be  signif 
cant  deformation  within  crustal  blocks  and  terranes 
[e.g.,  Davies,  1983). 

Several  models  have  been  proposed  for  the  high 
seismicity  in  northern  Alaska.  These  can  be  grouped 
into  those  involving  separate  tectonic  plates  [e.g.. 
Stone,  1983],  deformation  within  a  single  plate 
(Nakamura  et  al.,  1980;  Davies,  1983;  Biswas  et  al., 
1983 1,  or  a  combination  of  the  two  (this  study).  Stone 
( 19831,  among  others,  has  hypothesized  the  existence 
of  a  Bering-Arctic  plate,  which  includes  much  of 
Alaska  and  postulates  a  very  slow  moving  transform 
fault  along  the  Canadian  Arctic  Islands  coinciding 
with  a  diffuse  band  of  seismicity  extending  from  in¬ 
terior  Alaska  to  the  Beaufort  Sea  [Gedney  et  al., 
1977;  Biswas  and  Gedney,  1979], 

Focal  mechanisms  19  'this  study)  and  20  [Hase 
gawa  et  al.,  1979]  have  north-south  oriented  P  axes, 
perpendicular  to  the  trend  of  the  Nansen-Gakkel 
spreading  ridge,  and  nodal  planes  that  parallel  the 
continental  margin.  However,  the  low  level  of  seis 
micity  over  the  entire  margin  .nay  preclude  even  a 
very  slow  (less  than  1  cm/yr  [Stone,  1983])  deforma 
tional  margin. 

Nakamura  et  al.  (19801  show  stress  trajectories 
for  Alaska  based  on  the  direction  of  the  maximum 
horizontal  compression,  as  determined  from  the  ori 
entation  of  volcanic  dikes  and  active  faults  (Figure 
9).  Davies  ( 19>53I  introduced  the  concept  of  the  rigid 
indentor  as  applied  to  Alaska  to  explain  the  stress 
trajectories  (see  Molnar  and  Tapponmer  1 19751,  for 
an  explanation  of  the  concept  as  it  was  first  applied 
to  Indian  Asian  collision  tectonics).  Plasticity  theory 
predicts  a  fanning  of  the  stress  trajectories  from  the 
corners  of  a  rigid  indentor  being  pushed  into  a  plastic 
medium  For  Alaska  the  indentor  is  defined  as  the 
contact  area  between  the  Pacific  and  North  Amcri 


can  plates,  roughly  the  zone  between  the  trench  and 
the  50-km  contour  to  the  Wadati-Benioff  zone.  This 
model  explains  fairly  well  the  orientations  of  focal 
mechanism  pressure  axes  compiled  in  this  study. 
However,  the  amount  of  seismicity  in  northern 
Alaska  and  northwest  Canada  is  not  easily  explained 
in  terms  of  this  model  (Figure  9).  As  expected,  the 
most  active  seismicity  is  observed  along  the  main 
stress  trajectory  (west-nrrthwest  from  the  northern 
most  portion  of  the  indentor).  Yet  high  levels  of 
seismicity  are  also  observed  northeast  of  this  trend  in 
the  Yukon  Territory  with  little  seismicity  in  between 
(Figure  6). 

Central  to  any  satisfactory  hypothesis  of  the  pres 
ent  deformation  in  northern  Alaska  and  northwest 
Canada  must  be  an  explanation  of  the  prominence  of 
strike-slip  faulting.  Many  of  these  faults  show  right- 
lateral  offset  of  tens  to  hundreds  of  kilometers.  Table 
5  is  a  listing  of  the  major  right-lateral  strike-slip 
faults  with  their  corresponding  offsets,  if  known,  and 
the  ages  of  movement.  The  parallel  trends  of  the 
Tintina-Kaltag,  Denali,  and  Castle  Mountain-Lake 
Clark  fault  systems  arid  the  apparently  younger  mo 
tion  from  north  to  south,  together  imply  that  there 
has  been  a  consistency  of  process  through  time 
IGrantz,  1966;  Gedney  et  al.,  1974;  Gabrielse,  1985). 
The  Castle  Mountain- Lake  Clark  ( La  hr  et  al.,  19861, 
Denali  [Gedney  and  Estes,  1982],  and  Tintina  (this 
study)  faults  are  all  presently  showing  right  lateral 
strike-slip  motion.  This  suggests  that  the  Pacific- 
North  American  plate  boundary  is  distributed  across 
a  low-strain  rate  shear  zone,  approximately  500  km 
wide  and  extending  from  the  continental  margin  to 
at  least  the  Tintina  fault.  A  similar  suggestion  has 
been  made  by  Lahr  and  Plafker  (1980).  When  'be 
seismicity  and  focal  mechanisms  in  northern  AlasKa 
and  the  Yukon  are  seen  in  the  light  of  this  hypothe 
sis.  many  features  which  were  inconsistent  with  the 
other  models  are  explained.  The  seismicity  and  focal 
mechanisms  along  the  Eskimo  Lakes  fault  are  con 
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Fig.  9.  Model  proposed  for  deformation  of  continental  Alaska  and  northwest  Canada.  The 
exact  boundaries  of  the  "rigid  indentor”  (stippled)  are  speculative.  The  50-km  Wadati-Benioff 
zone  contour  is  from  Stephens  et  al.  [1984]  and  Stone  [1983].  Pacific-North  America  motion 
(RM-1)  and  convergence  rate  are  from  Minster  et  al.  [1974],  Stress  trajectories  (heavy  dashed 
lines)  are  from  Nakamura  et  al.  [1980],  selected  earthquakes  (M>6)  are  from  NEIS.  Small 
arrows  represent  the  P  and  T  axes;  the  large  arrows  represent  the  distributed  shear  zone  dis¬ 
cussed  in  the  text. 


sistent  with  compression  from  the  south.  The  seis¬ 
micity  occurs  at  a  sharp  bend  in  the  Rocky  Mountain 
front,  the  bend  being  a  possible  stress  concentrator. 
The  focal  mechanisms  in  interior  Alaska,  being 
mostly  strike-slip,  are  consistent  with  a  shear  system 
extending  into  the  area  from  the  southeast.  The 
northern  extent  of  the  shear  system  has  been  shown 
to  be  as  far  north  as  68°N  where  east-west  striking 
strike-slip  faults  cut  across  Brooks  Range  structures 
]  Dillon.  19851. 

The  Huslia  earthquake  described  earlier  and  the 
extensional  tectonics  on  the  Seward  Peninsula 
[Turner  et  al.,  1981;  Biswas  et  al.,  1983;  Lockhart, 
1984]  are  enigmatic  features  in  a  region  dominated 
by  strike-slip  faulting.  The  Huslia  earthquake,  the 
iargest  event  in  northern  Alaska  in  historic  times. 


must  fit  into  any  reasonable  tectonic  model.  Normal 
faulting  can,  however,  be  consistent  within  a  strike- 
slip  regime  if  the  concept  of  wrench  faulting  and  for 
mation  of  pull-apart  basins  is  invoked. 

Cenozoic  extension  on  the  Seward  Peninsula 
must  be  accounted  for  in  any  tectonic  model  of  north 
ern  and  western  Alaska.  This  extension  has  been 
studied  extensively  using  geologic  and  geophysical 
methods.  The  region  is  characterized  by  north-south 
extension  on  east-west  striking  normal  faults.  Ex 
tensive  Cenozoic  alkali  and  tholeiitic  basalts  [Turner 
et  al.,  1981]  and  a  zone  of  seismically  active  normal 
faults  in  the  southern  Seward  Peninsula  [Biswas  et 
al.,  1983]  are  manifestations  of  the  extensional  defor 
mation.  Lockhart  [1984]  concluded  that  the  exten 
sion  is  not  caused  by  localized  rifting  but  by  regional 


'I  nBLE  5.  Offsets  on  Major  Strike-Slip  Faults  in  Alaska  and  Northwestern  Canada 


Fault 

Age  of  Motion 

Offset 

Reference 

Castle  Mountain  -  Lake  Clark 

Early  Cretaceous  - 
Present 

tens  of  km 

Grantz  [  1966] 

Denali,  east  of  Mount  McKinley 

Cretaceous  -  Present 

400  km 

Turner  et  al.  11974; 

Denali,  west  of  Mount  McKinley 

100-115  km 

Grantz  [1966] 

Tintina.  in  Yukon  and  Alaska 

Paleozoic  and  Mesozoic, 
most  in  Cretaceous 

450  km 

Dover [1985] 

Tintina.  northern  British  Columbia 

Middle  Jurassic  -  Early 
Cretaceous 

+  900  km 

Gabrieise  [  1 985] 

Kaltag 

Late  Cretaceous  - 
Present 

140  km 

Patton  and  Hoare 
[1968] 
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extension  reactivating  ancient  faults.  This  is  based 
on  an  extensive  gravity  study  of  the  area,  on  compar¬ 
ison  of  the  Seward  Peninsula  to  other  rifts  around 
the  world,  and  on  the  high  proportion  of  alkali  versus 
tholeiitic  basalts,  indicating  a  differentiated  and  not 
a  primitive  mantle  source  for  the  basalts. 

The  Seward  Peninsula  regional  extension,  if  not 
caused  by  continental  rifting,  could  have  a  source 
similar  to  that  proposed  for  the  Huslia  area  of  the 
Koyukuk  basin:  wrench  faulting  in  a  compressive 
strike-slip  tectonic  regime.  A  problem  with  this 
model  is  the  immense  size  of  the  extensional  area 
compared  with  other  puli-apart  basins  [Aydin  and 
Nur,  19821.  Another  possibility  is  a  back  arc  tension 
al  regime.  However,  the  dimensions  of  the  area  and 
the  distance  to  the  arc  are  longer  than  are  commonly 
associated  with  back  arc  extension. 

Conclusions 

A  simple  P  wave  velocity  structure  of  a  35-km- 
thick,  6.1  km/s  crust  over  a  7.9  km/s  upper  mantle 
was  developed  for  interior  and  northern  Alaska.  The 
crustal  velocity  value,  determined  using  different 
techniques  in  the  Dali  City  area  and  eastern  Brooks 
Range,  is  consistent  with  continental  crust. 

The  Alaskan  North  Slope  is  apparently  aseismic, 
and  earthquakes  larger  than  magnitude  3.0  are 
evenly  distributed  within  and  south  of  the  Brooks 
Range  and  follow  the  curve  of  the  range  to  the  north¬ 
east.  This  seismicity  occurs  on  major  tectonic  faults 
and  also  within  terranes  or  crustal  blocks.  The  cen¬ 
tral  Kaltag,  northwest  Tintina,  central  Kobuk, 
Eskimo  Lakes,  Porcupine,  Dali  City,  and  Rampart 
faults  are  presently  seismirally  active.  The  largest 
earthquakes  (greater  than  magnitude  6.0)  occur  to 
the  south  of  the  Brooks  Range  on  a  west-northwest 
trend  parallel  to  the  main  stress  trajectory  that 
would  result  from  the  convergence  of  the  Pacific  and 
North  American  plates.  In  the  interior  of  Alaska 
they  occur  either  on  tectonic  faults  or  within  crustal 
blocks  at  the  intersection  of  lineaments. 

Focal  mechanisms  determined  for  earthquakes  in 
central  interior  Alaska  indicate  strike-slip  faulting 
with  north-northwest  oriented  compressional  axes. 
These  parallel  the  predicted  stress  trajectories  of  the 
Pacific  plate  indenting  into  North  America.  The  con¬ 
cept  of  large-scale  shearing  as  a  mechanism  of  defor¬ 
mation,  in  addition  to  the  compressional  interaction 
of  the  Pacific  plate  with  southern  Alaska,  can  ac¬ 
count  for  the  preponderance  of  strike-slip  mecha¬ 
nisms  on  major  crustal  faults.  We  propose  that  the 
Pacific-North  American  transform  boundary  in¬ 
cludes  a  zone  of  low  strain  rate  deformation  about 
500  km  wide  which  exhibits  right-lateral  strike-slip 
faulting,  'dost  of  the  relative  motion  is  on  the  Fair- 
weather  fault  (greater  than  5  cm/yr),  with  less  on  the 
Denali  fault  (about  1.0  mm/yr  [Horner,  1983]),  and 
probably  much  less  motion  on  the  Tintina  fault. 
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ADDITIONAL  EVIDENCE  FOR  DOWN-DIP  TENSION  IN  THE  PACIFIC 
PLATE  BENEATH  CENTRAL  ALASKA 

By  Larry  Gedney  and  John  N.  Davies 
Abstract 

Small  intermediate-depth  earthquakes  cluster  in  a  volume  about  20  km  in 
diameter  beneath  Gold  King,  a  landing  strip  on  the  northern  flanks  of  the  Alaska 
Range.  These  events,  almost  600  km  from  the  Aleutian  trench,  are  the  northern¬ 
most  which  can  be  identified  as  being  associated  with  the  Benioff  zone  extending 
inland  from  Cook  Inlet  west  of  Anchorage.  First-motion  observations  of  the  Gold 
King  events  at  local  stations  can  be  combined  to  determine  a  tensional  axis 
dipping  steeply  to  the  north  at  an  azimuth  just  west  of  north.  This  composite 
mechanism  suggests  that  the  slab  is  subjected  to  down-dip  tensional  forces, 
probably  as  a  result  of  gravitational  sinking. 

Introduction 

Gutenberg  and  Richter  (1941,  1954)  observed  that  the  volcanoes  associated  with 
the  Aleutian  Islands  persisted  northeastward  into  south-central  Alaska  accompa¬ 
nied  by  earthquakes  100  km  or  more  in  depth.  Benioff  (1954)  suggested  that  the 
focal  zone  of  the  intermediate-depth  events  described  an  inclined  slab  dipping  at 
28°  beneath  the  volcanic  arc.  Tobin  and  sykes  (1966)  and  Tarr  (1970)  showed  that 
the  intermediate-depth  earthquakes  extended  well  to  the  north  beyond  the  volcanoes 
to  the  vicinity  of  Mt.  McKinley  and  slightly  north  of  the  Alaska  Range.  Tarr 
(unpublished  data,  1971)  has  observed  that  these  events  were  smoothly  connected 
as  an  inclined  plane,  or  Wadati-Benioff  zone,  to  the  shallow  earthquakes  near  the 
Aleutian  trench.  Seismicity  maps  by  Meyers  (1976)  illustrate  this  relationship. 

During  the  past  20  yr,  regional  seismograph  networks  operated  by  the  University 
of  Alaska  and  the  U.S.  Geological  Survey  have  refined  the  picture,  showing  that  the 
inclined  seismic  zone  is  about  20  km  thick,  and  that  there  are  abrupt  changes  in 
trend.  One  of  the  changes  in  strike  occurs  near  Cape  Douglas  in  lower  Cook  Inlet, 
and  another  is  found  just  south  of  Mt.  McKinley  (VanWormer  et  al„  1974;  Davies, 
1975;  Lahr,  1975;  Estes,  1978;  Agnew,  1980). 

It  is  the  purpose  of  this  short  paper  to  describe  a  cluster  of  intermediate -depth 
earthquakes  just  south  of  Fairbanks  which  may  define  the  absolute  northern  tip  of 
the  subducting  north  Pacific  plate  in  North  America. 

The  Gold  King  Earthquakes 

We  refer  to  these  earthquakes  as  the  Gold  King  cluster,  after  a  nearby  landing 
strip  on  the  northern  flanks  of  the  Alaska  Range  (Figures  1  and  2).  Even  though 
this  grouping  is  quite  pronounced,  it  escaped  notice  for  a  number  of  years  because 
the  earthquakes  were  of  fairly  small  magnitude  and  well  separated  in  time.  Further, 
although  small  shallow  earthquakes  are  fairly  common  in  the  area,  it  had  not  been 
noticed  that  events  of  the  present  group  all  fell  within  a  narrowly  confined  inter¬ 
mediate-depth  range. 

The  cluster  lies  near  the  reported  epicenters  of  a  magntiude  65  earthquake  on  21 
January  1929  and  a  magnitude  7.0  event  on  16  October  1947.  The  epicenter  of  the 
first  event  was  about  64.0°N,  148. 0°W  (Davis  and  Echols,  1960),  although  this 
location  was  based  largely  on  felt  reports.  There  is  considerable  scatter  in  the 
various  locations  given  for  the  1947  earthquake  (St.  Amand,  1948;  Hodgson  and 
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Milne,  1951;  Davis  and  Echols.  1960;  Jordan  era/..  1968).  H.  Pulpan  (1986.  personal 
communication)  relocated  the  epicenter  to  64.35°N,  148.20° W,  but  there  is  no 
indication  given  by  any  of  these  workers  that  the  1947  earthquake  was  of  greater 
than  normal  depth,  and  it  seems  unlikely  that  the  Gold  King  earthquakes,  which 
are  mostly  100  km  or  more  in  depth,  can  be  associated  with  that  event. 


156°  154"  152'  150°  148*  146* 


154“  152"  150"  148'  146'  144'  142' 


Fig.  1.  Three-letter  station  designators  show  the  locations  of  the  seismographic  stations  utilized  in 
this  study.  The  two  upper  hemisphere  focal  mechanisms  on  the  left  were  obtained  by  Bhattacharva  and 
Biswas  1 1 979).  and  the  one  on  the  right  is  from  this  report.  Cross-section  A-A '  is  100  km  thick,  and  B- 
B  '  is  50  km  thick. 

Since  the  first  earthquakes  which  can  be  identified  as  belonging  to  the  Gold  King 
cluster  were  recorded  in  1969,  the  seismographic  networks  in  the  state  have  had 
many  different  configurations.  The  stations  which  have  operated  the  longest,  and 
the  ones  used  for  most  of  the  data  presented  in  this  report,  are  those  shown  in 
Figure  1. 
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The  better  solutions  for  the  Gold  King  earthquakes  are  those  available  from  1978 
through  1984,  and  86  of  these  are  plotted  in  Figures  2  and  3.  (In  addition,  Figure  3 
shows  all  regional  events  meeting  the  criteria  given  below.)  In  general,  the  solutions 
shown  in  Figure  2  fall  within  a  depth  range  of  100  to  130  km,  and  most  of  them  lay 


1  4  *•  \  4  T" 


8  p _ S  to _ <3  30 _ ?3  MH.ES 

S _ 0  a  IQ  ta  ?0  ?5  KILOMETERS 

Flu.  Epicentral  plot  of  86  earthquakes  of  the  Gold  King  cluster  located  between  1978  and  1984. 
WRH.  GKC,  and  MOK  are  local  seismographic  stations. 


within  a  radius  of  10  km  of  64.06°N,  148.H°W  (the  pronounced  center  of  the  group 
and  the  coordinates  of  an  event  described  later).  Although  some  vertical  scattering 
is  noted,  20  km  is  probably  a  good  approximation  for  the  dimensions  of  the  cluster 
as  a  whole,  because  depth  control  is  less  reliable  than  horizontal  control.  Magnitudes 
of  the  earthquakes  were  all  small;  the  largest  was  only  M/.  =  3.1. 
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Only  earthquakes  which  were  well-recorded  by  at  least  five  stations  have  been 
plotted  on  either  of  these  figures.  S-wave  arrivals,  which  were  typically  very  distinct, 
were  used  in  the  solutions  when  available.  The  crust  and  upper  mantle  structure 
given  by  Huang  and  Biswas  11983)  for  central  interior  Alaska  was  used  in  a  version 
of  the  U.S.  Geological  Survey’s  HYPOELLIPSE  location  program  (Lahr,  1984). 

The  events  plotted  in  Figures  2  and  3  were  required  to  have  rms  solution  errors 
of  less  than  0.5  sec.  Selected  events  used  for  a  composite  focal  mechanism  solution 
were  required  to  have  very  clear  first  arrivals  and  rms  solution  errors  of  less  than 
0.2  sec.  The  parameters  of  the  37  earthquakes  in  the  latter  group  are  given  in  Table 
1.  Nearly  all  fell  within  a  few  kilometers  of  the  area  of  densest  hypocentral 
concentration. 

Figure  1  shows  the  locations  of  the  two  northeast- looking  cross-sections  plotted 
as  Figure  3.  Section  A-A  ',  just  to  the  northeast  of  Mt.  McKinley,  clearly  shows  the 
characteristic  profile  of  a  Benioff  zone.  Just  beyond  this  to  the  northeast,  section 
L  .1 '  demonstrates  that  the  subduction  zone  has  all  but  died  out,  save  for  the 
remnant  that  the  Gold  King  cluster  represents. 

Focal  Mechanism 

Table  1  presents  a  summary  of  the  first  motions  that  are  plotted  in  the  focal 
mechanism  shown  as  Figure  4.  The  arrivals  at  most  stations  were  very  consistent 
over  the  years,  although  some  reversals  were  noted  at  stations  falling  near  a  nodel 
plane  (such  as  RDS  and  TOA). 

In  plotting  Figure  4,  the  azimuths  and  distances  given  to  the  individual  stations 
are  those  which  would  have  resulted  from  an  earthquake  occurring  at  coordinates 
64.06°N,  148.1 1°W  and  at  106  km  focal  depth.  This  location  falls  near  the  center 
of  the  cluster,  and  the  coordinates  are  those  of  the  earthquake  on  23  May  1978  at 
17:48  UTC. 

The  solution  parameters  are  given  in  Table  2.  From  Figures  3  and  4,  it  is  seen 
that  the  tensional  axis  lies  subparallel  to  the  surface  of  the  Pacific  lithospheric  slab 
at  an  azimuth  mid-wav  between  the  direction  of  plate  motion  (Minster  and  Jordan, 
1978)  and  the  plunge  vector  of  the  downgoing  slab  (Agnew,  1980).  Bhattacharya 
and  Biswas  (1979)  present  two  focal  mechanism  solutions  of  earthquakes  in  the 
adjacent  area  to  the  southwest  (on  5  March  1973  and  on  24  April  1974;  Figure  1)  in 
the  same  depth  range  (106  and  121  km,  respectively)  that  also  show  steeply  dipping 
T  axes,  although  the  orientations  of  the  null  axes  vary  considerably. 

Discussion 

To  the  south  of  the  present  study  area,  the  down-dip  direction  for  the  tensional 
axis  for  intermediate-depth  earthquakes  has  been  observed  within  the  Pacific  plate 
on  both  the  Alaska  Peninsula  (Akasche,  1968;  House  and  Jacob,  1983)  and  in  the 
Cook  Inlet  area  (Lahr,  1975;  Bhattacharya  and  Biswas,  1979;  Engle,  1982). 

However.  Enitle  (1982)  also  notes  that  there  is  an  alignment  of  the  P  axes  of 
shallow  and  intermediate  earthquakes  in  a  general  north-south  direction  around 
Cook  Inlet.  Further,  Pulpan  and  Frohlich  (1985)  argue  that  the  north-south  ori¬ 
entation  of  the  P  axes  in  the  region  between  59°N  and  60.5"N  is  stronger  than  the 
down-dip  alignment  of  the  T  axes.  They  conclude,  therefore,  that  the  focal  mecha¬ 
nism  data  indicate  a  “generally  north-south  horizontal  compression  regime"  beneath 
Cook  Inlet.  They  propose  that  the  transition  between  a  tensional  stress  regime 
south  of  59°  on  the  Alaska  Peninsula  to  the  compressional  regime  to  the  north  can 
be  attributed  to  segmentation  of  the  Pacific  plate. 
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This  segmentation  can  be  delimited  by  the  abrupt  30°  to  40°  bends  in  the  strike 
of  the  Benioff  zone  at  59°  N  and  63°  N.  A  consistent  interpretation  of  the  Gold  King 
data  would  be  that,  while  the  stress  regime  is  predominantly  compressional  in  the 
central  offset  section  betweeen  59°N  and  63°N,  down-dip  tension  is  the  dominant 
mechanism  in  the  slab  to  the  north  and  south  of  this  segment. 


64.06°N  148.1 1*W  H-106KM 

Fit;  4  Upper  hemisphere,  equal-area  plot  of  composite  focal  mechanism  solution  obtained  for  Gold 
King  cluster.  Solid  circles  and  open  circles  represent  initial  compressions  and  dilatations,  respectively. 
First  motions  at  individual  stations  are  tabulated  in  Table  1.  Down-dip  vector  is  taken  from  Agnew 
1 19801.  and  slip  vector  is  trom  Minster  and  Jordan  1 1978). 

TABLE  2 

Sommary  of  Composite  Focal  Mechanism  Solution  Parameters 

/>  Axi8  7*  Axis  H  Axis  Nodal  Plane  A  Nodal  Plane  H 

Azimuth  Plunge  Azimuth  Plunge  Azimuth  Plunge  Strike  Dip  -Strike _  Dip 

100  12  349  89  193  19  N26"E  o6°F.  NlO'W  40“  W 
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Summary . 


An  array  of  seismometers  was  deployed  in  the  area  of  the 
Yukon  Flats,  Alaska,  to  record  explosions  being  detonated 
along  the  Trans-Alaska  Pipeline  corridor.  The  purpose  of 
this  experiment  is  to  determine  the  crustal  structure  under 
the  Yukon  Flats  area,  and  under  interior  Alaska  by  combinina 
the  data  from  the  the  temporary  stations  with  data  from  the 
University  of  Alaska  permanent  seismometer  network.  The 
data  are  currently  being  analysed. 


Introduction 

The  original  proposal  was  designed  to  study  the  crustal 
structure  of  north-central  and  interior  Alaska  using  the 
Dali  City  and  Wood  River  sequences  of  earthquakes.  During 
the  course  of  this  work  it  was  discovered  that  the  US 
Geological  Survey  planned  to  detonate  a  number  of  large  (up 
to  6000  lb)  chemical  explosions  as  part  of  the  Trans  Alaska 
Lithosphere  Investigation  (TALI)  refraction  survey  across 
Alaska.  To  take  advantage  of  these  "artificial  earthquakes" 
as  an  aid  to  determining  the  crustal  structure,  the  proposed 
work  plan  was  modified  to  include  the  installation  of  a 
network  of  temporary  stations  in  the  Yukon  Flats  area.  The 
Yukon  Flats  area  was  selected  because  it  offers  the 
possibility  of  tying  in  the  existing  data  from  the  Dali  City 
event  (Estabrook  et  al.,  1987a)  with  data  obtained  by  the 
study  of  earthquakes  north  and  south  of  the  area  (Estabrook, 
1985,  Estabrook  et  al.,  1987b). 


Yukon  Flats  Experiment 

During  July  and  August  of  1987,  a  network  of  temporary 
seismometers  were  deployed  along  the  Yukon  and  Porcupine 
rivers  to  the  east  of  the  Trans-Alaska  pipeline  bridge, 
along  the  Steese  highway  to  Circle  and  at  the  villages  of 
Beaver,  Chalkyitsik  and  Birch  Creek  (figures  1  &  2).  Seven 
of  these  stations  were  three  component,  five  day  recorders 
loaned  to  us  by  the  USGS.  Five  of  these  were  deployed  along 
the  Yukon  river,  one  along  the  Porcupine  river  and  one  at 
Birch  Creek.  With  the  exception  of  the  Birch  Creek  station, 
they  were  all  deployed  by  river  boat.  The  remainder  of  the 
stations  were  a  mixture  of  smoked-drum  recorders  and 
seismometer  systems  designed  to  be  connected  directly  to  a 
telephone  line.  In  the  latter  case  the  seismic  signals  were 
transmitted  directly  to  recorders  in  the  Geophysical 
Institute.  Of  the  temporary  stations  two  were  manned,  and 
only  connected  during  the  scheduled  times  of  the  explosions. 
The  remaining  stations  were  fitted  with  telephone  answering 


devices.  These  were  set  up  so  that  after  ten  rings  the 
recorder  would  answer,  announce  itself,  switch  the  caller 
through  to  the  seismometer  signal,  switch  on  a  30  minute 
timer  then  switch  itself  out  of  the  circuit.  This  device 
saves  the  cost  of  long  long-distance  calls,  and  only 
requires  visits  to  install  and  dismantle  the  system. 

Because  the  Yukon  flats  area  consists  largely  of 
geologically  recent  river  deposits,  and  because  it  is  not 
possible  to  move  very  far  from  the  river  bank  for  those 
stations  installed  by  river  boat,  noise  surveys  were 
conducted  to  determine  whether  there  were  differences  in 
background  noise  between  different  types  of  river  bank.  It 
was  found  that  as  long  as  the  river  was  not  shallow  and 
turbulent,  distance  from  the  waters  edge  had  very  little 
effect  on  the  background  noise  levels.  However,  the  noise 
seemed  to  be  amplified  at  locations  were  the  trees  were 
small  or  non-existent  in  comparison  with  recently  cut  banks 
with  older  larger  trees.  We  ascribe  this  to  the  greater 
compaction  and  mechanical  rigidity  of  the  older  wooded 
soils. 

Out  of  the  13  temporary  stations  we  installed  we  obtained 
recordings  of  at  least  some  shots  from  all  but  2  of  them. 

In  these  two  cases,  equipment  failures  caused  the  loss  of 
record.  The  data  from  the  five-day  recorder  stations  are 
currently  being  played  back  by  the  USGS  at  Menlo  Park.  The 
stations  recorded  on  smoked  drums  or  via  telephone  lines  to 
the  Geophysical  Institute  have  been  read,  and  time-distance 
plots  are  being  constructed. 

In  addition  to  the  temporary  network,  the  fixed  network 
stations  run  by  the  Geophysical  Institute  Seismology 
Laboratory  (figure  1)  were  also  used  to  record  the  shots. 

All  of  the  shot  signals  being  received  at  the  Geophysical 
Institute  were  recorded  on  both  film  recorders  and  directly 
on  the  dedicated  MASSCOMP  computer  system.  Tables  1  and  2 
give  the  locations  of  the  explosions  and  the  numbers  of 
stations  in  the  fixed  network  which  recorded  them.  To  date, 
only  first  arrival  signals  have  been  studied,  but  it  is 
planned  to  use  the  digital  recordings  to  construct  composite 
record  sections.  The  data  from  the  USGS  recorders  will  also 
be  digitised  from  the  tapes  and  included  in  the  composite 
sections . 


1987  Dali  City  Event 

During  the  first  sequence  of  USGS  shots  there  was  a 
magnitiude  5.5  event  just  to  the  north  of  the  Trans-Alaska 
Pipeline  crossing  of  the  Yukon  River.  As  can  be  seen  from 
figure  1,  it  is  almost  on  line  with  the  stations  deployed 
along  the  Yukon  and  Porcupine  rivers.  This  event  was 
recorded  at  four  of  the  three  component  five-day  recorder 


sites  (two  other  sites  having  fallen  to  either  battery 
problems  or  damage  by  bears).  Since  these  data  can  all  be 
digitised,  this  should  provide  an  excellent  data  set  for 
coda  and  attenuation  studies,  and  provide  additional 
information  on  the  crustal  structure,  though  it  is  beyond 
the  scope  of  this  project  to  pursue  these  studies  very  far. 


Results  to  date. 

As  mentioned  above,  the  data  from  the  explosions  and  the 
Dali  City  event  recorded  during  the  summer  of  1987  are 
currently  being  processed.  In  addition  to  this  data,  models 
of  posible  crustal  velocity  structures  based  on  the  very 
limited  data  available  in  the  literature,  and  a  re-anlvsis 
of  explosion  data  from  the  1985  USGS  refraction  line,  have 
been  assembled  (Stone  e"  al.,  1987).  We  plan  to  have  the 
preliminary  analysis  of  the  data  completed  by  the  American 
Geophysical  Union  meeting  to  be  held  in  San  Francisco  in 
December  1987. 


Future  Explosion  Studies 

The  USGS  plan  to  continue  their  refraction  experiments  in 
1988,  with  lines  along  the  Fairbanks  -  Circle  road  (Steese 
highway)  ,  and  the  Livengood  -  Manley  road.  Vie  are  planning 
to  deploy  the  five-day  recorders  again  to  obtain  ray-paths 
that  cross  the  Yukon  Flats  from  the  south  and  east  to  the 
north,  to  complement  the  dominantly  south  to  north  and  east 
to  west  paths  obtained  during  the  1987  experiment. 

Starting  in  October  1987,  the  US  Army  intend  to  detonate  a 
series  of  large  (9000  lb)  chemical  explosions  on  the  south 
flank  of  the  Alaska  Range  near  Donnelly  Dome.  These  energy 
sources  can  be  utilised  to  investigate  the  crustal  structure 
north  of  the  Alaska  Range  to  Fairbanks.  We  plan  to 
supplement  the  fixed  seismometer  network  with  about  ten 
dial-up  stations.  A  tentative  distribution  of  stations  is 
shown  in  figure  2.  This  distribution  should  allow 
constraints  to  be  placed  on  the  crustal  structure  beneath 
the  Tanana  Valley  and  beneath  the  Alaska  Range. 
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Figure  1. 

A  map  showing  the  combined  temporary  (cir 
fixed  network  (stars)  seismometer  station 
1987  Yukon  Flats  and  interior  Alaska  crus 
experiment.  Also  shown  are  the  shotpoint 
USGS  refraction  lines,  and  the  Dali  City 
1987  (diamond) . 


cled  sta 
s  used 
tal  stud 
s  (  +  )  f  o 
event  of 


rs)  and 
for  the 
les 

r  the  three 
24  July 


Figure  2. 

A  map  showing  the  distribution  of  seismometer  stations, 
larger  shot  points  and  the  major  faults  of  Alaska.  The  shot 
points  are  labelled  by  shotpoint  number  followed  by  the  size 
of  the  explosion  in  thousands  of  pounds.  The  large  circles 
are  of  100  and  200  km  radius  about  the  proposed  winter 
explosions.  The  small  open  triangles  are  tentative 
locations  for  dial-up  systems. 
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Table  1. 

Shotpoint  locations  and  times  of  detonation 


SHOT 

SHOT 

DATE 

SHOT  TIME 

LAT 

LONG 

SIZE 

NO 

PT 

(DAY : YR) 

(HR:MIN : SEC) 

(H) 

(VO 

(LBS) 

1 

19 

205:87 

08:00:00.016 

61 

01.83 

145 

15.68 

6000 

2 

44 

205:  87 

08:02:00.012 

63 

38.96 

145 

53.59 

1000 

3 

49 

205  :  87 

08:04:00.007 

65 

17.81 

146 

25.16 

6000 

4 

3 

205:  87 

08:06:00.010 

62 

54.65 

145 

29.75 

3000 

5 

8 

205:87 

08:08:00.008 

61  . 

,56.02 

145 

17.54 

6000 

6 

42 

205:87 

08:10:00.005 

63 

14.59 

145 

38.52 

500 

7 

48 

205:87 

10:00:00.007 

64 

54 .39 

146 

21.76 

4000 

8 

1 

205:87 

10:02:00.005 

63 

20.78 

145 

42.13 

2000 

9 

45 

205:  87 

10:04:00.012 

63 

47.35 

145 

50.74 

2000 

10 

41 

205:87 

10:06:00.010 

63 

02.28 

145 

30.82 

1000 

11 

6 

205:87 

12:00:00.076 

62 

29.51 

145 

28.07 

4000 

12 

2 

205  :  87 

12:02:00.010 

63 

09.87 

145 

32.00 

2000 

13 

46 

205:87 

12:06:00.012 

63 

55.33 

145 

47.67 

3000 

14 

43 

205:  87 

12:07:30.005 

63 

29.38 

145 

50.90 

1000 

15 

54 

206:87 

10:00:00.008 

64 

22.03 

146 

07.74 

3000 

16 

46 

206:87 

10:02:00.010 

63 

55.33 

145 

47.67 

3000 

17 

52 

229:87 

08:02:00.014 

63 

52.14 

145 

12.01 

3000 

18 

51 

229:87 

08:04:00.008 

63 

32.32 

144 

00.38 

4000 

19 

50 

229:87 

08:04:00.009 

62 

52.80 

143 

05.99 

6000 

20 

56 

229:87 

08:06:00.012 

64 

32.96 

146 

47.95 

2000 

21 

66 

229:87 

08:08:00.008 

65 

41.69 

149 

11.32 

6000 

22 

55 

229:87 

08:10:00.014 

64 

28.84 

146 

35.51 

500 

23 

53 

229:87 

10:00:00.012 

64 

01.72 

145 

34.92 

2000 

24 

57 

229:87 

10:02:00.012 

64 

38.67 

147 

02.00 

3000 

25 

73 

229:87 

10:04:00.013 

64 

22.02 

146 

12.45 

700 

26 

59 

229:87 

12:00:00.009 

65 

04.72 

147 

41.77 

4000 

27 

47 

229:87 

12:02:00.010 

64 

13.86 

146 

00.81 

2000 

28 

50 

232:87 

08:00:00.016 

62 

52.80 

143 

05.99 

6000 

29 

69 

237:87 

08:00:00.007 

66 

18.29 

150 

25.84 

4000 

30 

62 

237:87 

08:02:00.006 

65 

19.85 

148 

18.44 

2000 

31 

54 

237:87 

08:04:00.007 

64 

22.03 

146 

07.74 

5000 

32 

57 

237:87 

08:06:00.157 

64 

38.67 

147 

02.00 

4000 

33 

65 

237:87 

08:10:00.007 

65 

34.47 

148 

56 .31 

4000 

34 

67 

237:87 

10:00:00.012 

65 

47 . 36 

149 

24.89 

3000 

35 

61 

237:87 

10:02:00.006 

65 

16.89 

148 

07.84 

1000 

36 

70 

237:87 

10:04:00.005 

66 

36 .12 

150 

59.19 

4000 

37 

64 

237:87 

10:06:00.011 

65 

28.76 

148 

41.70 

1000 

38 

66 

237:87 

12:00:00.015 

65 

41.69 

149 

11.32 

2000 

39 

60 

237:87 

12:02:00.005 

65 

11.93 

148 

04.76 

2000 

40 

59 

237:87 

12:04:00.241 

65 

04.71 

147 

41.79 

3000 

41 

63 

237:87 

12:06:00.011 

65 

24 . 35 

148 

25.14 

2000 

42 

74 

238:87 

08:00:00.005 

65 

03.42 

150 

11.26 

3000 

43 

63 

238:87 

12:02:00.010 

65 

24.35 

148 

25.14 

1100 

44* 

61 

240:87 

22:00:00.007 

65 

16.89 

148 

07.84 

45* 

59 

240:87 

22:02:00.009 

65 

04.71 

147 

41.79 

46* 

60 

240:87 

22:30:00.007 

65 

11.93 

148 

04.76 

47* 

60 

241:87 

01 : 30:00.007 

65 

11.93 

148 

04.76 

48* 

59 

241:87  01:32:00.009 

65  l 

04.71 

147  ■ 

41.79 

49* 

61 

241:87  02:00:00.007 

65  : 

16.89 

148  l 

07.84 

*  SHOT  ST7.F  T  S  1  D  0  -  ?  fl  D  r.BS 


Table  2. 

Sequential  shot  number,  shotpoint,  size  of  explosion  and  the 
number  of  fixed  network  seismometers  that  recorded  it. 


SHOT  NUMBER 

SHOT  POINT 

SIZE 

(LBS) 

STATIONS* 

1 

19 

6000 

12 

2 

44 

1000 

19 

3 

49 

6000 

12 

4 

3 

3000 

23 

5 

8 

6000 

12 

6 

42 

500 

4 

7 

48 

4000 

2 

e 

1 

2000 

9 

45 

2000 

1 

10 

41 

1000 

16 

11 

6 

4000 

8 

12 

2 

2000 

17 

13 

46 

3000 

2 

14 

43 

1000 

15 

54 

3000 

8 

16 

46 

3000 

17 

52 

3000 

15 

18 

51 

4000 

21 

19 

50 

6000 

•  - 

20 

56 

2000 

10 

21 

66 

6000 

7 

22 

55 

500 

11 

23 

53 

2000 

15 

24 

57 

3000 

17 

25 

73 

700 

11 

26 

59 

4000 

10 

27 

47 

2000 

16 

28 

50 

6000 

18 

29 

69 

4000 

5 

30 

62 

2000 

7 

31 

54 

5000 

22 

32 

57 

4000 

14 

33 

65 

4000 

7 

34 

67 

3000 

8 

35 

61 

1000 

-  - 

36 

70 

4000 

4 

37 

64 

1000 

7 

36 

66 

2000 

.  . 

39 

60 

2000 

5 

40 

59 

3000 

7 

41 

63 

2000 

1 

42 

74 

3000 

9 

43 

63 

1100 

3 

4  4  *  * 

61 

5 

4  5  *  * 

59 

5 

46** 

60 

47** 

60 

2 

48*  * 

59 

3 

49** 

61 

5 

*  NUMBER  OF  STATIONS  THE  BLASTS  WERE  OBSERVED  AT 
**  SHOT  SIZE  100-200  LBS 
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DETERMINATIONS  IN  THE  FAIRBANKS,  ALASKA  AREA 
BASED  ON  ARTIFICIAL  EXPLOSION  EXPERIMENTS 


APPENDIX  VII 


Presented  to  the  Faculty  of  the  University  of  Alaska 
in  Partial  Fulfillment  of  the  Requirements 
for  the  Degree  of 

MASTER  OF  SCIENCE 

Jean  M.  Johnson,  BJ. 

May  19W 


ABSTRACT 


Ten  explosions  which  occurred  in  the  Fairbanks  area  are  used  to  determine, 
first,  how  accurately  seismic  e rents  can  be  located  using  the  computer  program 
HYPOELLIPSE;  second,  the  factors  which  influence  how  well  the  events  are  lo¬ 
cated;  and,  third,  which  of  the  four  proposed  velocity  models  best  approximates 
the  Fairbanks  area.  The  average  hypocentral  error  is  4.9  km,  the  average  epicen- 
tral  error  is  2.6  km,  and  the  average  depth  error  is  3.2  km.  The  most  important 
factors  influencing  the  location  accuracy  are  station  configuration,  reading  errors, 
distance  to  closest  station,  and  the  velocity  model  used.  The  HYPOELLIPSE 
error  indicators  are  shown  to  be  smaller  than  the  actual  errors.  The  proposed 
velocity  models  are  all  faster  than  the  actual  structure  of  the  Fairbanks  area;  how¬ 
ever,  the  actual  errors  show  that  the  events  are  fairly  accurately  located.  The  data 
is  too  limited  to  propose  a  new  velocity  model. 


